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I.  INTRODUCTION 


In  this  paper  we  study  analytically  the  electromagnetic  radiation  field  produced  by  the  cloud  of 
accelerating,  non-relatlvistic  elections  induced  at  a  planar  photoelectron  emitting  surface  in  vacuum  by 
a  photon  pulse  of  arbitrary  time  profile  and  small  spatial  cross-sectional  area  incident  upon  that 
surface,  deriving  scaling  relations  for  the  radiated  fields  (hence  radiated  power,  energy,  and  spectral 
content)  explicit  in  the  pulse  and  surface  parameters. 

The  dynamics  of  die  electron  cloud  has  been  of  interest  in  connection  with  both 
electromagnetic  pulse  (EMP)  phenomena  [1-5]  and  laser  plasma  jets  [6-7]  and  has  been  studied  both 
analytically  »nd  numerically.  Those  studies  are  concerned  with  either  steady-state  dynamics  [1,3,6]  — 
applicable  when  the  election  plasma  period  is  much  shorter  than  die  characteristic  photon  pulse  width 
—  or  with  transient  dynamics  [2, 4,5,7],  applicable  when  the  above  condition  does  not  hold.  Those 
authors  who  obtain  complete,  strictly  analytical  results  for  the  transient  case  [2,7]  assume  1-D  planar 
geometry,  monoenergetic  emission  into  vacuum,  and  a  "no-screening"  approximation  [1 ,2,4,7].  Under 
these  assumptions,  they  solve  Maxwell's  equations  (simply  Gauss’s  law  in  this  situation)  for  the 
electric  field  in  die  space  charge  region,  coupled  to  Newton’s  law  for  election  motion  as  well  as  to  a 
continuity  equation  for  the  election  density,  to  obtain  a  self-consistent  solution  to  die  dynamics 
problem.  This  program  produces  formal  expressions  for  the  election  density  and  velocity  fields  as  a 
function  of  the  single  space  coordinate  and  time;  these  expressions  could  be  used  in  principle  to  derive 
similarly  formal  expressions  for  die  radiated  fields  of  interest  to  us  here,  although  this  has  not  been 
done  in  die  literature  as  far  as  we  know.  However,  from  our  point  of  view  there  is  a  fundamental 
difficulty  with  these  formal  expressions:  They  can  be  made  to  yield  expressions  to  die  radiated  fields 
explicit  in  the  pulse  and  surface  parameters  only  in  the  two  special  cases  of  photon  pulses  which  are 
either  constant  in  time  or  linearly  increasing  in  time.  Indeed,  these  two  special  cases  are  the  only  ones 
offered  in  [2]  and  [7]  as  illustrations  of  their  general  formalism.  The  specific  difficulty  is  this.  The 
treatment  of  die  particle  dynamics  is  done  from  a  Lagrangian  (in  the  fluid  sense)  point  of  view  while 
the  electron  density  and  velocity  fields  are  inherently  Eulerian  (as  is  usually  die  case  in  the 
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electrodynamics  of  extended  charge  distributions).  The  formal  expressions  for  the  density  and  velocity 
fields  require  a  translation  from  Eulerian  to  Lagrangian  coordinates  in  order  to  be  expressed  explicitly 
in  terms  of  the  pulse  and  surface  parameters  ([2].  Eqs.  (18)  and  (36);  [7],  Eqs.  (16M18))  and  this 
translation  can  only  be  done  exactly  analytically  in  the  two  special  cases  indicated  above. 

In  order  then  to  achieve  generality,  we  abandon  the  usual  technique  of  using  the  (Eulerian) 
charge  and  current  densities  to  compute  die  radiation  fields;  rather,  we  compute  these  fields  by  the 
novel  technique  of  directly  summing  contributions  to  the  fields  over  individual  electron  trajectories, 
i.e.,  we  do  "Lagrangian  (in  the  fluid  sense)  electrodynamics".  (We  point  out  that  this  technique  is 
novel  only  because  it  is  being  used  as  a  strictly  analytical  tool  -  the  technique  of  summing 
contributions  over  individual  particles  has  been  used  previously  in  electromagnetic  particle-in-cell 
(PIC)  codes.)  This  summing  process  results  in  integral  expressions  for  the  electric  and  magnetic  fields 
at  large  but  finite  (i.e.,  "finitely-remote")  distances.  True  radiation  quantities,  however,  are  obtained 
only  in  die  limit  of  die  field  point  going  to  infinity;  we  thus  demonstrate  the  existence  of  this  limit  and 
obtain  asymptotic  radiation  quantities  in  which  all  die  integrations  have  been  fully  carried  out  and 
which,  in  addition,  are  explicit  in  the  pulse  and  surface  parameters.  Because  we  wish  to  obtain  this 
limit  rigorously,  die  entire  paper  of  necessity  takes  on  a  somewhat  mathematical  flavor  in  order  to 
support  from  die  outset  this  major  goal.  We  have,  however,  relegated  most  of  the  mathematical  detail 
-  which  involves  only  elementary  techniques  -  to  four  appendicies  so  that  they  may  be  avoided,  for 
die  most  part,  by  those  readers  who  find  them  of  lesser  interest 

In  summary  then,  we  derive  expressions  for  die  asymptotic  radiation  fields  based  upon  a  2-D 
cylindrical  electron  cloud  of  small  radius  whose  dynamics  is  described  by  the  aforementioned  1-D 
planar  model  employing  the  monoenergetic  emission  and  no-screening  approximations.  Since  in  the 
1-D  model  die  magnetic  field  is  strictly  absent  then  that  model  yields  electrostatic  interactions 
between  the  electrons;  hence,  application  of  die  1-D  dynamical  model  to  our  2-D  electron  cloud  is 
tantamount  to  assuming  only  electrostatic  interactions  and  so  errs  in  neglecting  small  magnetic  field 
and  retardation  effects  in  the  dynamics.  Such  effects  are  not  neglected  in  our  radiation  treatment 
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however;  indeed,  we  require  a  small  spatial  cross-sectional  area  for  the  photon  pulse  so  that  we  may 
adequately  represent  the  retarded  time  in  die  radiation  integrals. 

In  Section  II  we  describe  the  photon  pulse  more  completely.  In  Section  III  we  discuss 
election  cloud  dynamics,  presenting  in  detail  only  those  features  unique  to  our  treatment.  In  Section 
IV  we  derive  expressions  for  the  finitely-remote  fields;  these  lead,  in  Section  V,  to  the  asymptotic 
fields  which  in  turn  lead,  in  Section  VI,  to  all  asymptotic  radiation  quantities  of  interest  Finally,  in 
Section  VII  we  present  illustrations  of  our  general  formulation  for  five  specific  pulses:  constant,  linear 
ramp,  triangular,  parabolic,  and  sin2.  Additionally,  Appendicies  A  and  B,  C,  and  D  contain  die 
mathematical  complements  to  Sections  HI,  IV,  and  V  respectively. 
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II.  PHOTON  PULSE  DESCRIPTION 

The  photon  pulse,  normally  incident  on  die  emitting  surface,  is  taken  to  be  one  of  circular 
spatial  cross-section  with  spot  radius  a,  and  photon  frequency  v.  This  restriction  to  circular  cross- 
section  is  not  necessary,  since  we  will  ignore  edge  effects  in  the  cylindrical  electron  cloud,  but  seems 
most  natural.  We  assume  that  photoelectrons  are  emitted  spatially  uniformly  over  the  extent  of  the 
spot,  each  with  one  and  the  same  non-relativistic  speed  v0  >  0  in  a  direction  normal  to  the  locally 
planar,  smooth  emitting  surface,  with  average  photoelectron  yield  (quantum  efficiency)  over  the  spot 
denoted  by  Y(v)  (electrons/photon).  (Here,  "locally"  means  over  a  region  comparable  to  that  of  the 
spot)  Monoenergetic,  normal  emission  of  electrons  is  certainly  not  the  case  physically  in  general 
where,  even  for  a  normally  incident  pulse  of  monochromatic  photons,  the  shape  and  width  of  the 
emitted  photoelectron  kinetic  energy  distribution  function  is  highly  dependent  upon,  and  varies  greatly 
with,  the  nature  and  condition  of  the  emitting  surface  as  well  as  with  the  incident  photon  frequency 
[8];  but,  as  we  have  previously  pointed  out,  monoenergetic  (normal)  emission  is  a  standard  assumption 
in  the  literature.  We  also  assume  that  the  incident  photon  intensity  is  not  high  enough  to  produce  a 
plasma  at  the  surface;  and  we  do  not  include  contributions  to  the  radiation  field  of  induced  currents  in 
case  die  emitting  surface  is  a  conductor  [6]. 

The  full  width  of  the  pulse  in  time  is  taken  as  n/ft,  where  ft  >  0  is  an  angular  frequency 
(sec1),  while  die  pulse  intensity  at  the  emitting  surface  is  given  by 

/  (t)  =  Af(t),  0£t£  n/ft,  (2.  l) 

where  A  >  0  is  an  amplitude  (photons/m2-sec  (we  use  MKSA  units  throughout])  and  f  is 
dimensionless,  with  0  £  f(t)  £  1,  f(t)  *  1  for  at  least  one  t,  and  f(t)  =  0  for  at  most  t  €  {0,  n/ft}.  It  is 
convenient  to  introduce  a  dimensionless  time  variable  given  by  s  =  ftt;  in  terms  of  s,  we  d*fi»** 
functions  /  and  g  by 
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/( s)  =  /(s /ft)  2  0,  g(s)  =  f(s/ft)  £  0, 


0  <,  s  £  re, 


(2.2) 


so  that 


/(s)  =  Ag(s), 


0  £  s  <i  n. 


(2.3) 


Hie  central  construct  of  our  model  turns  out  to  be  j  /(tOdt'  *  A  ff(t  7)dt 7  in  terms  of  which 

0  0 

ail  physical  qimntitw  of  interest  may  be  naturally  written.  It  will  often  be  more  convenient  to  use  a 
dimensionless  multiple  of  this  integral:  we  define 


a  k 

G(s)  *  J  g(s')d s'  /J  g(s/)ds/.  0  £  s  £  it. 


(2.4) 


so  that 


or  *  i  «D 

G(ftt)  «  J  g(s)ds  /  J  g(s)ds  *  J  f(t')dt'  /  J  f^t  ,)dt/.  OSt^Ji/Q. 


(2.5) 


Also,  we  define 


tf)  * 

N(ft)  *  J  f(t')dt'  «  (1/ft)  J  g(s)d s 


(2.6) 


so  that 


j  f(t')dt'  *  N(ft)G(ftt),  0  £  t  £  tt/ft; 
further,  we  exclude  as  trivial  the  zero  pulse  so  that  N(ft)  >  0.  Therefore 


(2.7) 
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0  £  G(S)  £  1  for  0  £  s  £  x  and  G(0)  *  0,  G(x)  -  1. 


(2.8) 


We  also  require  that  G''  exists  on  [0,  it]  (taking  one-sided  limits  at  0,  x);  to  ensure  this  we  demand 
that  g  be  continuous  on  [0,  x]  (or  f  on  [0,  x/fl]),  and  so  is  properly  integrable  there,  and  in  that  case 
we  have 

G'(s)  -  g(syON(a>  « /(syAON(fl)  k0,  0  £  s  <;  x.  (2.9) 

Hence  G  is  non-decreasing  on  [0,  x]  and,  being  continuous  there,  takes  on  every  value  between  0  and 
1  (by  Eq.  (2.8)).  In  fact,  since  g  vanishes  only  at  most  at  0  and  x  then  G/  is  strictly  positive  on 
(0,  x)  and  hence  G  is  strictly  increasing  on  [0,  x). 

We  point  out  that  we  could  easily  extend  our  results  to  improperly  integrable  pulses  g  which 
are  continuous  only  on  (0,  x)  but  for  which  lim  g(s)  *  «  or  lim  g(s)  *  «»  (or  both),  so  that  we  could 

•-KT  wr 

treat,  for  example,  g(s)  *  A  [(x-s)*“-X’w];  and  to  pulses  of  infinite  width  which  are  continuous  on 

(0,  «■»)»  so  that  we  could  treat,  for  example  g(s)  *  Ae  s  2:  0  (s  >  0) ;  but,  on  physical  grounds, 

we  see  little  reason  to  do  so.  We  could  also  allow  g  to  have  discontinuities  in  the  interior  of  its 

interval  of  definition  so  that  we  could  treat,  for  example,  g(s)  =  A- 11^’  ^  i  <  ^  .  This  last 

l  1,  It  TIL  S  S  S  X 

generalization  only  makes  mathematical  arguments  more  cumbersome  and  provides  essentially  no 
additional  physical  insight  We  choose  not  to  pursue  any  of  these  possibilities  here. 
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ID.  ELECTRON  CLOUD  MOTION 

A.  Kinematics 

At  any  instant  t  £  0,  we  think  of  the  electron  cloud  as  a  circular  cylinder  of  radius  a 
comprised  of  a  continuum  of  infinitesimally  thin,  planar,  parallel  electron  disks  (or  "sheets"),  each 
having  been  emitted  at  some  instant  t0  €  [0,  t]  from,  and  each  remaining  forever  parallel  to,  die  planar 
emitting  surface,  with  precisely  one  such  sheet  having  been  emitted  for  each  t0  e  [0,  t].  It  follows 
that  each  sheet  may  be  unambiguously  labeled  by  the  time,  tD,  at  its  emission.  Observe  that  the  total 
charge  contained  in  die  sheet  emitted  during  infinitesimal  time  interval  8tc  at  tD  is  (forever) 

8q(g  -  -eY (y)KxjKa36ta  -  -eY(v)A0N(0)G  '(Ot^na  38t0, 

where  e  is  the  magnitude  of  the  electron  charge,  and  that  this  sheet  has  infinitesimal  thickness 
5z  =  Vofito  so  that  also 

8q(tJ  *  -eYi(tJn(fbz/vo.  (3.2) 

Further,  if  we  let  a  z-axis  coincide  with  the  symmetry  axis  of  the  cylinder  and  choose  as  positive  the 
direction  of  electron  emission,  then  we  denote  die  position  at  time  t  of  the  sheet  labeled  by  t0  as 
Z(t0;  t)  for  (tc;  t)  e  [0,  n/O)  x  (0,  »);  clearly 

Z(t„;  t)  £  0  for  all  to,  t,  (3.3) 

die  emitting  surface  being  located  at  z  =  0.  Figure  1  depicts  the  electron  cloud  and  the  emitting  surface 
at  any  particular  instant,  L 

We  assume  that  the  magnitude  of  the  residual  positive  charge  on  the  emitting  surface  is 
precisely  the  same  as  the  magnitude  of  die  totality  of  negative  charge  in  the  electron  cloud;  that  is,  we 
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CHARGE  CLOUD 


EMITTING  CHARGE 

SURFACE  SHEET  t. 


0  Z(t^  t)  Z«(t) 


Figure  1.  Election  cloud  description  for  dynamics. 


8 


assume  that  neutralization  of  this  residual  positive  charge  by  the  flow  towards  it  of  negative  charges 
originating  in  die  emitting  surface  -•  in  case  it  is  metallic  ~  is  negligible  during  time  intervals  of 
interest  here.  Further,  since  the  thickness  of  the  positively  charged  layer  is  on  the  order  of  10*  m  [8], 
we  may  safely  assume  that  it  forms  a  surface  charge  distribution. 

B.  Dynamics 

In  general,  the  net  electromagnetic  force  at  time  t  on  the  charge  sheet  labeled  by  ^  is 

F(to; t)  -  l&KgtafizlSz  JJ  {E(r.  0. Z(t0; t),  t)  ♦  V(t„;  t)  x  B(r,  d,  Z(t0;  t),  t)}dS  (3.4) 

(we  assume  any  one  cylindrical  coordinate  system  whose  positive  z-axis  is  the  one  previously  specified 
following  Eq.  (3.1))  where  V(ta;  t)  is  the  Instantaneous  sheet  velocity,  E(r,  d,  Z(te;  t),  t)  and 
B(r,  A,  Z(to;  t),  t)  are  the  instantaneous  electric  and  magnetic  fields  at  the  area  element  dS  located  at 
(r,  0)  on  the  sheet,  and  n  =  no2.  Equation  (3.4)  holds  since  the  sheet  has  uniform  charge  density 

PL(t„)  -  5q(to)/n5z  -  (l/nv0)(8q(toy8t0].  (3.5) 

As  alluded  to  earlier,  we  approximate  the  foil  electromagnetic  field  E  +  (V  x  B)  by  the  net 
electrostatic  field,  E**,  resulting  from  the  charge  in  the  cloud  not  on  sheet  to  plus  that  in  the  residual 
surface  charge  distribution.  This  net  electrostatic  field  is  given  by 

Z-<» 

E  "(r,  d,  Z(t0;  t),  t)  *  -k(l/eo)  J  dC p(C,  t)  ,  (3.6) 

Z(V«) 

where  k  is  the  unit  vector  in  the  +z  direction,  Z»„(t),  1 2  0,  is  the  z-coordinate  of  the  leading  edge  of 
the  charge  cloud,  and  p(C,  t)  £  0,  (£,  t)  e  [0,  is  the  z-dependent  volume  charge  density 
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characterizing  the  charge  cloud;  hence  we  have 

ZJl) 

F(t0;  t)  -  -*{5q(toy£j  J  d{p(t,  t)  -  fiqOJ  •  E  “(r,  0,  Z(tc;  t),  t).  (3.7) 

Z^i«) 

The  reader  should  note  that  because  die  z-component  of  die  electrostatic  force  at  any  surface  element 
dS  located  on  charge  sheet  te,  due  to  all  the  other  charge  in  the  system  not  on  sheet  ts,  is  independent 
of  r  (the  radial  coordinate  on  die  sheet),  and  because,  in  addition,  all  die  charge  on  sheet  t0,  other  than 
that  belonging  to  dS,  exerts  no  z-directed  force  on  dS,  then  the  z-component  of  die  force  on  dS  due  to 
all  other  charges  in  the  system  is  independent  of  r;  hence,  sheet  tD  remains  forever  strictly  planar  in 
this  approximation.  On  the  other  hand,  while  all  the  charges  not  on  sheet  t0  exert  no  net  radial  force 
on  dS  in  our  approximation  (Eq.  (3.7)),  die  charges  on  sheet  tc  and  not  belonging  to  dS  do  exert  a  net 
radial  force  on  dS  and.  after  some  time,  significant  radial  expansion  will  occur.  We  assume  negligible 
radial  expansion  during  the  time  intervals  being  considered  here. 

Hie  equation  of  motion  of  the  cenler-of-mass  of  charge  sheet  to  is 

6nKt0)[(a1Z/dtJ)(t0;t)]  -  F(t„;t)  (3.8) 

where  8m(to)  is  the  mass  of  the  sheet  and  F(t0;  t)  given  by  Eq.  (3.7).  That  is, 

Z*.(0 

(yz/atXa)  *  /  <*Cp(C.t)  *  o  .  0.9) 

« f.;0 

Of  course,  this  equation  is  valid  only  for  those  times  t  for  which  the  sheet  t0  has  not  previously 
returned  to  die  emitting  surface,  namely  for  t  satisfying 

t0  S  t  <  !♦,  (3.10) 
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where  t*  is  given  by 


t*  «  sup  (T  e  (t0,  -)  |  Z(tc;  t)  >  0  for  ail  t  e  (t„;  T)}  (3.11) 

(possibly  t*  ■  «-);  note  t*  >  t0  by  Eq.  (3.12)  below.  In  other  words,  given  t0  €  [0,  n/Q],  to  solve 
the  equation  of  motion  we  must  find  a  function  t  Z(t0;  t)  and  a  t*  e  (to; «]  such  that 
(1)  Z(tc;  t)  satisfies  Eq.  (3.9)  for  t  e  [t0,  t*)  and  initial  conditions 

Z(te;  t0)  =  0  and  0Z/atXto;  ta)  =  vc  >  0 ;  (3.12) 


and 

(2)  t*  and  Z(t0;  t)  satisfy  Eq.  (3.11). 

Outside  tiie  time  interval  specified  by  Eq.  (3.10)  the  following  must  hold: 

Z(t0;  t) «  0  *  (dZ/3tXt„;  t)  =  (d’Z/dt’X t„;  t)  if  te  (— ,  t„)  u  (tj,  ~).  (3.13) 

(We  will  specify  values  for  these  functions  at  (tQ;  t*),  in  case  t*  <  «>,  later  on.)  In  particular,  it 
follows  from  this  and  Eq.  (3.12)  that  t  **  (3Z/dtXt„;  t)  is  discontinuous  at  t  ■  t0  so  that  (d2Z/dt2Xtp;  t) 
does  not  exist  there;  hence  the  second  derivative  in  Eq.  (3.9)  should  be  interpreted  to  be  the 
appropriate  one-sided  quantity  at  t  =  to.  The  physical  interpretation  of  Eq.  (3.9)  is  that  the 
instantaneous  acceleration  of  sheet  tD  at  time  t  depends  only  upon  the  total  amount  of  charge  in  the 
cloud  between  the  sheet  itself  and  the  cloud’s  leading  edge,  being  proportional  to  that  charge. 

Equation  (3.9),  as  an  equation  for  unknown  Z(t0;  t),  is  not  readily  solvable  since  Z(t0;  t) 
occurs  as  a  limit  of  integration  and  since  Z^tt)  and  p((,  t)  are  also  unknown.  In  order  to  proceed  it 
is  useful  to  make  one  additional  simplifying  assumption,  namely,  that  charge  sheets  do  not  pass 
through  one  another.  This  is  the  so-called  "no-screening"  or  "no-charge-sheet-crossing"  (NCSC) 
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approximation  If  such  is  the  case,  then  the  integral  in  Eq.  (3.9)  is  independent  of  time  and  is  given 
by  (using  Eqs.  (32),  (2.1),  and  (2.7)) 

a  f  dC(XC,t)  -  |8q(t0)  -  -eYaJfctjdt. 

0  .  (3.14) 

-  -eY Aaf  f(t0)dt0  -  -eYAflN(Q)G(nt0)  (NCSC). 

O 

Hence  Eq.  (3.9)  becomes 

02Zfit2Xto;t)  -  ^esAYN(Qym,e0]G(Qt0)  (NCSQ,  (3.15) 

valid  for  t„  £  t  <  t*  (with  t*  yet  to  be  determined).  This  is  the  required  equation  of  motion  of  charge 
sheet  t0.  It  implies  that  a  given  sheet  t0  experiences  a  constant  (in  time)  negative  acceleration  and, 
furthermore,  if  tj  >  to  then  sheet  tj  experiences  a  more  negative  accderatkm  than  does  sheet  tc 
(since  G  is  strictly  increasing  on  [0.  x]). 

We  may  now  integrate  Eq.  (3.15)  twice,  subject  to  initial  conditions  given  by  Eq.  (3.12),  to  get 
0Z/at)(to;t)  -  v„  -  [e 2AYN(Q)hnce0]G(Qto)(t  -t„)  (3.16) 


and 


Z(tc;t)  -  vo(t-tc)  -  (l/2)[e2AYN(nym,E0]G(Gt0Xt  -t,,)2. 


(3.17) 


Defining 

■  eHAY/voym,e0 , 


(3.18) 
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where  (Up  is  the  plasma  frequency  corresponding  to  maximal  electron  emission  density  AY/v0,  we  have 


(dZ/dt)(t0;t)  -  vJl-^N(Q)CX£*0)(t-to)) 


(3.19) 


and 


Z(t0;t)  -  v0(t-t„)[l  - (l/2)o^ N(Q)G(£lt0)(t - 10>]  ; 


(3.20) 


also 


(0f2Z/3t*)(to;  t)  -  -v0a£N(Q)G(ftto)  .  (3.21) 

We  reiterate  dial  these  equations  are  valid  only  for  t„  £  t  <  t*  and  for  NCSC. 

To  find  t*.  we  observe  that  file  motion  of  charge  sheet  t„  is  clearly  that  of  an  object  rising  and 
falling  under  "gravitational”  acceleration  of  magnitude  voa^N(fl)G(flto)  and,  as  such,  Z(to;  t) 
increases  monotonically  from  z  *  0  at  t  =  t„  to 

*  -  Z*(U  ■  d/2)vo  [l/<ajN(«)G(nt0)]  -  (1/4)^770(00  (3.22) 

at 


tw(0  ■  tQ  ♦  (l/(^N(Q)G(Qt0)]  -  t„  ♦  (1/2)770(00  >  tc 
and  then  decreases  monotonically  to  z  =  0  at 

WO  ■  t„  ♦  T/G(Gt0)  >  t0 


(3.23) 


(3.24) 
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where 


T  m  2/«>}N(Q)  >  0  (3.25) 

is  the  round-trip  duration  for  die  sheet  launched  at  t0  *  xin  (assuming  no  t0  cutoff  —  see  the  next 
subsection).  It  follows  that  t*  given  by 

to  -  WU  (126) 

satisfies  Eq.  (3.11),  with  t*  *  - if  tc  =  0,  so  Eqs.  (3.19)  -  (3.21)  are  valid,  given  tc  €  [0,  ic/fl],  for 

te  [t0,t0>77G(flto)).  (3.27) 

In  summary,  we  may  combine  Eqs.  (3.19)  -  (3.21)  with  Eq.  (3.13)  to  arrive  at  expressions  for 
position,  velocity,  and  acceleration  of  the  charge  sheet  with  label  to  e  [0,  n/i 2],  for  all  t  e  [0,  •»),  as 
follows: 


Z(t#;t) 


vc(t  -gn-  r'cxnga  -t„)j,  if  t  e  i(t„) 

0,  ift€l-(t0). 


(3.28) 


vjl  -27~‘G(Gt0)(t  -t0)),  if  t  €  I(tD) 

0,  if  t  e  I  "(t0)  , 


and 


(3.29) 
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A(t0;  I) 


-2(v0/7)G(ft0),  if  t  e  I(t0) 

0.  if  t  e  I'(t0)  , 


(3.30) 


where  we  have  used 

KO  ■  lO  t0  +  77G(Qt0)]  -  [t0,t*J 


(3.31) 


and 

1(0-  (-.«-)  \I(t„).  (3.32) 

We  will  also  denote,  for  future  use, 

I°(0  -  <*„•  +  7TCKOO)  -  K  O  .  (3.33) 

Note  that  Z(ta;  • )  is  continuous  at  t„  and  t*  but  neither  V(ts;  • )  nor  A(tD;  • )  is  at  either,  further,  we 
have  arbitrarily  elected  to  define  these  latter  two  maps  at  t*  so  that  they  are  left-continuous  there  (we 
have  already  arranged  that  they  be  right-confinnous  at  te).  Also, 

V(t0;t)  -  (3Z/9t)(to;t)  and  A(t0;t)  -  (aPZ/3t*)(t0;t)  ,  t  *  to,C  (3.34) 

It  follows  from  above  and  the  previous  paragraph  that 

|V(t„;t)|  ta  v0,  teKO  (3.35) 

so  that  the  motion  is  always  non-relativistic  if  it  is  initially  so. 
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C.  The  Cutoff 

The  results  expressed  in  Eqs.  (3.28)  -  (3.30)  are  constrained  by  the  NCSC  requirement  Since 
thu  constraint  was  iBipn«wi  on  die  already  existing  equation  of  motion,  Eq.  (3.9),  which  does  not  a 
priori  exdude  charge  sheet  crossing,  then  we  must  insure  that  our  complete  model,  including  NCSC, 

Is  self-consistent  That  such  self-consistency  does  not  automatically  obtain  in  general  may  be  seen  as 
follows.  First  define  quantity  i.  by 

T„  ■  TJ7,  G]  ■  sup  fr  e  (0,  x/O]  |(dtra/dt0)(T)  £  0  for  all  x  e  (O.T]}.  (3.36) 

where  t^to)  Is  given  by  Eq.  (3.24);  note  that 

(<Vdt0)(t0)  -  1  -  12nG/(Gt0yGJ(Qt0))  (3.37) 

exists  on  [0,  x/Q]  since  G'  does,  that  t,  £  x/Q,  and  that  tD  satisfies 

G/(CI0)/GJ(£if0)  -  1 IQT  (3.38) 

whenever  ^  <  x/fl  Now  consider,  in  Rg.  2,  two  typical  pictures  in  the  to  - 1  plane  of  die  set 

U  HU  *  1CU1  of  arguments  for  which  sheet  motion  is,  according  to  Eqs.  (3.28)  -  (3.30),  non- 

stationary.  The  behavior  of  t^  near  tB  «  0  follows  from  limt^fU  “  00  which  in  turn  follows  from 

limGffltJ  *  qo,  while  the  behavior  dsewhere  follows  from  die  fact  that  G  is  strictly  increasing. 
v-*0* 

Then  it  is  clear  that  NCSC  does  not  hold  in  situations  represented  by  Rg.  2b. 

To  rectify  this  situation,  we  have  two  options:  We  may  either  restrict  ourselves  to  considering 
only  pulses  for  which  die  behavior  in  Rg.  2a  is  representative,  which  behavior  is  characterized  by 
LHS  (3.37)  £  0  for  tc  e  [0,  x/Q];  or  we  may  consider,  in  addition,  pulses  represented  by  Rg.  2b  if  we 
cut  off  from  consideration  in  the  dynamics  (and  hence  subsequently  in  the  electrodynamics)  those 
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Eqs.  (3.28H3.30). 


-  t.  +  77G(OU  =  UO 


t  =  t. 


77G(OtJ  >  0 


<  x/Q 


charge  sheets  with  to  >  to.  (In  fact,  we  could  cut  off  somewhat  beyond  to  [2]  but  choose  not  to  do  so 
since  we  are  interested  in  all  t  >  0.)  Elimination  of  these  late-time  charge  sheets  does  not  at  all 
influence  the  dynamics  of  the  stained  earlier-time  sheets  with  labels  t*,  ^  t0  since  die  dynamics  of  any 
particular  charge  sheet  is  influenced  only  by  those  charge  sheas  emitted  earlier  than  it  (in  the  NCSC 
approximation).  Further,  when  ^  does  not  differ  greatly  from  x/fl  then  using  this  cutoff  is  not  an 
unreasonable  approximation.  In  the  literature,  this  issue  is  addressed  only  in  [7]  where  the  first  opdon 
is  chosen.  We  choose  die  second  since  it  is  more  general  than  the  first,  including  the  latter  as  a 
special  case.  Naturally,  we  call  ^  the  cutoff  on  to.  Finally,  since  we  have  based  our  strategy  on  Fig. 

2,  we  should  point  out  that  the  two  possibilities  illustrated  there  are  not  only  typical  but  in  fact  general 
except  for  the  possible  occurence  of  points  of  inflection  or  intervals  of  constancy  in  t^to),  which 
features  do  not  affect  the  above  argument  Because  such  features  will,  however,  have  an  influence  on 
our  radiation  considerations,  we  demonstrate  in  Appendix  A  die  existence  of  a  class  of  pulses  which 
have  die  interval-of-constancy  feature. 

While  it  is  clear  that  elimination  of  charge  sheets  with  to  >  4  is  a  necessary  condition  for  self¬ 
-consistency  (when  we  wish  to  consider  all  t  >  0),  it  is  not  a  priori  dear  that  it  is  also  suffident  That 
is,  we  may  have  to  further  limit  those  to’s  which  we  admit  to  the  charge  cloud  in  order  to  insure 
consistency  with  the  NCSC  constraint.  We  now  settle  this  issue,  showing  that  this  condition  is  in  fact 
also  sufficient.  To  that  end,  we  let  label  to  be  admissible  if  charge  sheet  to,  during  its  entire  flight, 
never  spatially  coinddes  with  any  other  earlier-launched  charge  shed;  i.e.: 

Definition  3.1:  ta  e  [0,  Tt/fl]  is  admissible  if  for  every  t/  e  [0,  te)  and  for  every  t  e  P(t0)  n 
r(t0')  we  have  Z(t^ ;t)*Z(t0;  t). 

It  is  clear  from  this  definition  that  tc  =  0  is  always  admissible  (vacuously);  physically  this  may 
also  be  seen  since  Z(0;  t)  =  vet  for  any  pulse.  In  fact,  as  we  demonstrate  in  Appendix  B,  the 
following  is  true: 


18 


Theorem  3.2 :  (i)  t,,  >  0 

(it)  If  t0  6  [0,  Tj  then  t0  is  admissible;  further,  [0,  tj  is  the  largest  possible 
interval  of  admissible  t0’s  which  contains  to  =  0. 

In  summary,  the  NCSC  requirement  has  led  us  to  limit  (when  tD  <  n/Ci)  those  charge  sheets, 
labeled  by  t0,  which  we  admit  to  die  charge  cloud.  As  a  consequence,  if  we  define 

D[T,G]«  U  [<t0}  x  I(t0)]  (3.39) 

«.«  I0.T.1 

as  well  as  denote  its  complement  and  interior,  respectively,  by 

D-[T,G]  -  {[O.Tj  x  (-,-)}  \  D[T, G]  (3.40) 

and 

D°[T,G]  -  (J  [{tj  x  I°(t0)j  ,  (3.41) 

».«  (O.T.) 

then  the  maps  (t0,  t)  Z(tc;  t),  V(t0;  t),  A(t0;  t)  of  Eqs.  (3.28)  -  (3.30)  will  henceforth  be  considered 
to  have  only  domain  [0,  Tj  x  (—»,  <*>)  with  their  values  on  D[T,  G]  being  given  by  the  respective  first 
lines  of  Eqs.  (3.28M3.30)  and  their  values  on  D~[r,  G]  being  0.  In  particular,  they  will  no  longer  be 
considered  to  be  defined  for  values  of  t„  in  (?Q,  n/ft).  We  will  present,  in  Section  VII  examples  of 
pulses  with  T0  =  n/Ci  and  of  pulses  with  T„  <  nJQ.. 

D.  Charge  and  Current  Densities 

Since,  as  pointed  out  in  the  Introduction,  our  formulation  of  electrodynamics  is  driven  by  the 
inability  to  translate  analytically  from  the  Eulerian  to  the  Lagrangian  description  of  die  electron  cloud 
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charge  density,  p(z,  t),  and  current  density, 


J(z,  t)  =  p(z,  t)v(z,  t)  (3.42) 

(where  v(z,  t)  represents  the  charge  cloud  velocity  at  (z,  t)),  we  devote  some  discussion  to  this  issue. 
Consider  first  fire  charge  density.  Given  z  >  0  and  t  >  0,  we  expect  that 

p(z,t)  -  0ro/az)(z,  t)pL(x0(z,  t))  (3.43) 

where  t0  *  t0(z,  t)  is  "the"  label  that  satisfies 

Z(t0;t)*z  (3.44) 


for  pL  (tD)  given  by  Eq.  (3.5)  and  Z(to;  t)  given  by  Eq.  (3.28).  This  expectation  is,  of  course, 
basically  correct  but  we  must  be  careful  to  ensure  that  the  prescription  given  by  Eq.  (3.44)  is  well- 
defined,  i.e.,  that  such  a  t0,  given  by  a  differentiable  (wit  z)  tQ(z,  t),  exists  and  is  unique.  Indeed,  Eq. 
(3.44)  has,  in  general,  many  solutions  for  to  (given  z  and  t)  and  no  valid  physical  solution  for  t0 
whenever  z  >  vQt  (since  such  solutions  would  be  acausal). 

It  is  easy  to  settle  the  existence  and  uniqueness  issue  for  Eq.  (3.44).  To  this  end,  let  D,[T,  G] 
be  the  set  defined  in  die  z-t  plane  -  see  Fig.  3  -  by 

D,[T.  G]  ■  {(£,  x)|  there  exists  t0  e  [0,  T„]  such  that  Z(tc;  t)  *  £}.  (3.45) 

Qearly  D,[T,  G]  *  0.  Then  the  following,  whose  proof  is  provided  in  Appendix  B,  is  true: 
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Figure  3.  The  extendon  of  Fig.  2b  to  three  dimensions  to  show  some  representative  points  on  the 
surface  Z<t0;  t),  (t0;  t)  €  D[r,  G];  is  a  typical  value  of  t0. 


2 


Theorem  3.3:  Let  z,  t  >  0.  If  (z,  t)  €  D,(7,  G]  then  the  equation  Z(tc;  t)  =  z  has  a  unique 
solution  in  the  interval  [0,  Tj;  and  if  (z,  t)  e  Dj[T,  G]  then  p(z,  t)  =  0.  Fiirther,  if  z  >  0  and  t  £  0 
then  p(z,  t)  *  0. 

The  existence  of  dxjdz  at  points  (z,  t)  e  D,[T,  G]  follows  from  the  Implicit  Function  Theorem;  but, 
since  we  will  not  base  our  electrodynamics  on  p(z,t)  and  J(z,t),  we  will  not  give  die  details  here. 
(However,  the  essentials  of  the  argument  are  contained  in  a  portion  of  the  proof  of  Theorem  C.4.) 

Having  found  die  appropriate  t„  =  x„(z,  t)  required  for  charge  density  in  Eq.  (3.43),  we  then 
find 


v(z  ,t)  »  kV(to(z,  t);  t) ,  (3.46) 

where  V(t0;  t)  is  given  by  Eq.  (3.29),  and  finally  find  J(z,  t)  from  Eq.  (3.42)  so  that 

J(z,  t)  «  k(dxo/9zKz,  t)pL(xo(z,  t))V(T0(z,  t);  t) .  (3.47) 


In  practice  one  finds  t0  by  solving  Eq.  (3.44)  directly,  either  analytically  or  -  necessarily,  in 
many  cases  -  numerically  and  by  then  judiciously  selecting  the  correct  tc  from  among  the  in  general 
many  roots  (all  but  one  of  which,  however,  lie  outside  [0,  TJ).  However,  p  and  J  are  not  generally 
explicitly  analytically  available  via  Eqs.  (3.43)  and  (3.47)  because  Eq.  (3.44)  for  t„  is  not  generally 
solvable  analytically.  This  state-of-affairs,  then,  drives  our  choice  of  technique  for  doing,  in  the  next 
section,  charge  cloud  radiation  electrodynamics. 
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IV.  RNTTELY-REMOTE  FIELDS 

A.  Background 

The  usual  methods  for  studying  radiation  theory  for  macroscopic  charge  and  current 
distributions  employ  a  field-theoretic  (Eulerian)  description  of  the  distributions.  In  this  approach,  one 
assumes  that  p(x,  t)  and  J(x,  t)  are  available  and  from  them  calculates  the  retarded  potentials  at  space¬ 
time  points  (x,  t)  of  interest  according  to 

♦H(x,t)  -  (l/4jce0)JJJ dV(x')p(x',  t  -(|x  -x'|/c))|x-x/|-‘  (4.1) 

and 

Ab(x, t)  -  (1/4xeocj)JJJ dV(x')J(x',t-(|x -x'  |/c))|x  -xT  (4.2) 

where  the  integrals  extend  over  any  volume  containing  all  the  sources  (x  and  jf  are  position  vectors 
with  respect  to  some  arbitrary  origin);  die  fields  F.(x,  t)  and  B(x,  t)  then  follow  from 

E(x,t)  *  -V<Kx,  t)  -  (9A/dt)(x,  t)  and  B(x,t)  *  V  x  A(x,t).  (4.3) 

All  electrodynamic  quantities  of  interest  follow  from  E  and  B;  in  particular,  radiation  quantities  follow 
from  those  parts  of  E  and  B  falling  off  with  distance  as  l/|x|  as  opposed  to  l/|x|2. 

As  pointed  out  in  detail  at  the  end  of  the  previous  section,  the  above  Eulerian  description  is 
not  adequate  and  we  must  resort  to  another  approach.  This  other  approach  is  to  use  a  particle- 
theoretic  (Lagrangian)  description  of  the  macroscopic  charge  and  current  distributions  which  give  rise 
to  the  radiation  fields.  This  method,  while  common  for  studying  the  fields  of  a  single  moving  charge, 
is  not  generally  used  for  analytical  treatments  of  macroscopic,  continuously  distributed  aggregates  of 
charged  particles;  but  it  is  the  method  we  use  here.  This  Lagrangian  point-of-view  will  allow  us  to 
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obtain  analytical  expressions  for  all  the  usual  radiation  quantities  of  interest  without  sacrificing  the 
generality  we  require.  In  this  approach,  a  macroscopic  charge  distribution  is  considered  to  be  a  union 
of  "infinitesimar  charge  elements,  5q(Q.  with  label  I  in  some  index  set,  St,  where  5 q(D  has 
position,  velocity,  and  acceleration  denoted  by  X(t,  t),  X(t,  t),  and  Xfl;  t)  respectively.  (The  *  and 
■  should  be  regarded  here  only  as  distinguishing  notation  -  we  do  not  mean  to  imply  that  t  «•  X(t;  t) 
is  twice  differentiable  on  all  of  (-»,  «*).)  The  potentials  corresponding  to  Eqs.  (4.1)  and  (4.2)  are  then 
given  by 

♦L(x, t)  «  (1/4xeo) f  dq(f)l*  - Xfl;  t '(t  x.  t))|-‘[l  -0(1; x, t)]“*  (4.4) 


and 


Al(x,  t)  -  (l/4xE0c2)j  dq(f)X(t,  1 7(l;  x,  t))  j  x  -  Xfl;  t  '(I;  x,  t))  |  '*[  1  -0(1;  x.t)]'1 


(4.5) 


where 


0(1;  x,  t)  -  (l/c)X((;t7fl;x,t))  -R[x,X(t;t'(«;x,t)))  ,  (4.6) 

R[x,X((;t/((;x,  t))]  -  [x-Xfl;t7fl;x,t))]  /(x-Xfl;t7fl;x,t))  |  (4.7) 

is  the  unit  vector  pointing  from  Xfl;  t'fl,  x;  t))  to  x,  and  t'fl;  x,  t)  is  the  retarded  time  which  is  the 
solution  of 


t  -t 7  -  (1/c)  |x  -Xfl;  t #)  |.  t 7  *  t,  (4.8) 

given  (,  x  and  t  The  Adds  E  and  B  are  then  once  again  determined  by  Eq.  (4.3).  (It  is  assumed  that 
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Eq.  (4.8)  has  a  unique,  well-behaved  solution  for  t',  given  t,  x,  t;  we  will  have  more  to  say  about  this 
later.)  In  case  the  distribution  is  a  single  point  charge  q  (so  that  Sf  =  {1}  and  the  label  f  may  be 
suppressed,  and  |&Kl)h  4)  ***  potentials  and  are  the  Lidnard-Wiechert  potentials  of  the  point 
charge. 

Upon  performing  foe  operations  indicated  in  Eq.  (4.3)  for  E  and  B  we  find 

E(x,  t)  -  (1/4*6,,)  f  dq(f)  [£.(1;  x,  t  '(&  x,  t))  ♦  £.(*  x,  t  '(I;  x.  t))]  (4.9) 


and 


B(x,t)  «  (l/4xe0c)  f  dq(f){R[x, X(f; t y(f; x, t))l  x  [EJit, x, t '{t, x, t))  ♦  £,(l;x,t/(l;x.t))]} 


i 


(4.10) 


where 


«,(l;x.t'(tx)t))  -  [R -(l/c)X(f; t ')]  /  Y*(i 1 0[  1  -0(1; x, t)]J |x  - X(t t ')  |\  (4.11) 

-  (r  x  {[R-(l/c)X(l;t')l  x  Xftt')})  /  [1  -0(t;x,t)]3|x  -X(t;t') 1. 

(4.12) 

Kit')  -  (1  -  [|X(l;r/)|/c]J}'l/l;  (4.13) 


in  die  above  we  have  abbreviated  by  ft  die  quantity  which  is  represented  more  fully  in  Eq.  (4.7),  and 
by  tf  the  retarded  time  t'ft  x,  t).  The  subscripts  "v"  and  "a"  on  £  stand  for  "velocity”  and 
"acceleration";  the  velocity  fields  fall  off  as  l/|x  -  X|:  and  are  essentially  static  in  character,  while 


the  acceleration  Adds  fall  off  as  l/|x  -  X|  and  are  the  dynamic  radiation  Adds.  In  what  follows  we 
will  be  interested  only  in  die  radiation  Adds 

E,(x.t)  -  (1/4 «„)  f  dq(l)E,(t;  x.  t  '(tx,  t))  (4.14) 


and 


B.(x,t)  -  (l/4xeoc)  J dq(D{ R{x, X(t; t '(t x. t»]  x  ff.ttx.t'ttM))}. 


(4.15) 


The  Foyndng  vector  associated  with  this  electromagnetic  field  is 


S.CM)  «  (e#cJ)E„(x,t)  x  B,(x,t)  (wattahn1)  (4.16) 

so  did  die  instantaneous  radiation  energy  flux  area  density  at  x  in  die  direction  i  is  S.(x,  t)  -f  and 
the  instantaneous  radiation  energy  flux  angular  density  at  x,  |x|»0,  in  the  direction  f  (also  called 
instantaneous  radiated  power  per  unit  solid  angle)  is 

(8P/84*)(x,t)  -  S,(x,t)  -\x\H  -  cocat  •|*|E1(x,t)  x  |x|B.(x,t)  (watls/sr)  (4.17) 

(where  solid  angle  is  denoted  by  Y).  The  total  radiated  energy  angular  density  at  x  in  die  direction  t 
is 

(5W/8YXX)  -  J(8P/5Y)(x,t)dt  -  J(5P/8Y)(x,t)dt  (Joules/sr)  (4.18) 

0  •m 

and  die  spectral  intensity  of  this  radiation  is  given,  with  CD  £  0,  by 
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(52I/8*F&o)(x,a>)  -  2£0cJ|x|*Re[i  -E^(*.o))  x  B*(x,<o)]  (joules/sr-hz) 


(4.19) 


where 


cats  -  «r-/- 


and  *  denotes  complex  conjugate;  Eq.  (4.19)  follows  from  the  requirement 


(4.20) 


(8W/8Y)(x)  «  f  dto(5JI/54/8<o)(x,  co).  (4.21) 

0 

B.  The  Fldd  integrals 

Up  to  this  point,  everything  we  have  done  in  this  section  has  been  very  general.  We  now 
return  to  the  particulars  of  our  problem  which  are  illustrated  in  Fig.  4.  In  that  figure,  each 
infjnitKrimai  charge  element  of  foe  electron  cloud  is  labeled  by  the  label  ta  of  foe  shed  on  which  it 
(forever)  lies  and  by  its  plane  polar  coordinates,  (r,  0),  on  that  sheet  which  also  remain  constant  in 
time,  (r  remains  constant  because  of  our  assumption  of  negligible  radial  expansion  while  0  remains 
constant  because  of  fire  cylindrical  symmetry  of  the  charge  cloud.)  Thus,  our  set  of  labels  for  the 
charge  dements  is 


a  -  {(Vr,0)  |t„  e  [0.TJ,  r  e  [0,a],  0  e  l-*,it)  if  r  *  0,  and  0  -  0  if  r  -  0).  (4.22) 


Using  Eqs.  (3.1)  and  (3.18),  we  have 


8q(to,r,0)  -  [8q(toyiw2]rdrd0  -  ^YAGN(n)G/(tt(>)5tordrd0 
-  Km.e(/e)voa^QN(a)G/(Qto)fitordrd0  ; 


(4.23) 


Figure  4.  Ejection  cloud  description  for  electrodynamics.  The  symbol  t  may  be  interpreted  as  present 
or  retarded  time,  as  needed.  The  vector  x  is  taken  to  lie  in  the  first  quadrant  of  die  x-z 
plane  and  outside  die  cloud  (bence  x  •  1  >  a). 
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XO^r,  fr;t)  -  i  r  cosfr  ♦  J  r  sinfr  ♦  k  Z(t„;t), 


(4.24) 


XO^r.fct)  -  k  V(t„;t) 


(4.25) 


and 


Xd^r.ftt)  -  k  A(t0;  t) 


(4.26) 


where  Z,  V  and  A  are  given  by  Eqs.  (3.28H3.30)  with  t  e  (—,  «>). 

We  now  write  out  expressions  for  the  fields  E,(x,  t)  and  B,(x,  t)  based  on  Eqs.  (4.14)  and 
(4.15).  Since  we  have  assumed  that  v0  Is  non-relativistic,  it  follows  from  Eq.  (3.35)  that  so  is  V(te;  t) 
for  all  (tD,  t)  e  [0,  TJ  x  (—», «);  herre  re  take  in  Eq.  (4.12) 


[1  -  0(to,r.fr;x.t)r3  -  1  ♦  3(l/c)Xa0.r.d;t/(to.r.fr;x.t))  -R[x,  X<t0.r.fr;t'(t6,r.fr;x,t))]. 


Denoting 


Xd(^o»  0 


1, 

0, 


if  M)  e  Dtr.G] 
if  (t^t)  e  DtT.G] 


(427) 


(4.28) 


and  noting  from  Eqs.  (4.25)  and  (4.26)  that  X(t„  r,  d;  t)  x  X(t„,  r.  fi;  t)  *  0  we  find,  for  t  £  0, 


S  9 

^(x,t)»K J dt„G '(QtJGiQtJ Jdrr  Jddx0(t0;t/){[(R*k)R-k][l  +  W(t0;t')]  /Ix-X^r.fct')!} 


(4.29) 
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where  t7 » t/(t0,  r,  0;  z,  t)  is  retarded  time  (possibly  t7  <  0)  satisfying 
t  -  t7  -  (l/c)Jz  -  X(t„,r,fl;t7)|  .  t7  £  t; 


(4.30) 


also 


ft  -  ft[x,X(t0,r,  0;t7)], 

(4.31) 

W(tB;t7)  -  3(l/c)V(t0;  1 7)(ft  •  k) 

(4.32) 

and 

K  ■  (m./nec  2)fl(vo/772 . 

(4.33) 

Likewise,  for  t  £  0, 

T.  •  * 

B.(x,t)«-<K/c)  f dt0G7(OtoX3(«to) fdrrf<H>xD(t0;t'){(&  x  k)[l 

0  0  -t 

+W(t0;t7)]  /|x-X(t0,r,  fi;t7)|}. 

(4.34) 

Note  that  die  factor  t„;  t7)  ensures  that  sheet  t0  contributes  to  the  fields  only  if  it  is  in  flight  at  the 
retarded  time  t 7  i.e.,  if  t7  e  I  (tc).  It  is  dear  at  this  point  that 

E,(x,t)  ■  0  ■  B.(x,t)  for  0  £  t  £  |x|/c 

(4.35) 

although  we  will  re-derive  this  result  in  a  more  systematic  way  later. 

To  proceed  further,  we  concentrate  upon  the  integrands  in  die  above  expressions.  We  remark 
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that  by  symmetry  we  may,  without  loss  of  generality,  take  x  to  lie  in  the  first  quadrant  of  the  x-z 
plane  so  that 


x*j  «  0.  (4.36) 

Also,  the  observation  point  x  is  assumed  to  lie  outside  the  cloud;  hence 

4  *  x*I  >  a  >  0.  (4.37) 

With  a  view  towards  addressing  the  components  of  E,  and  B,  we  compute,  using 

C  ■  x-k,  (4.38) 

the  components  of  the  integrands  to  be  (abbreviating  X  =  X(tc,  r,  d;  t7)) 

[(ft-k)fc-i]  -  (k-i)  *  (l/|x-X|J)  [(x-X)-k]  [(x-X) •!] 

-  (l/|x— X  |*)  [^-C(X-i)  -  $(X*k)  ♦  (X'i)(X*k)] 

(4.39) 

*  (l/|x-X|2)  [4C~Cr  cost)  -  §Z(t0;t)  +  Z(t0;t)r  cosd] 

-  (l/|x-X|2)  (Z(tc;  t)  -  CKr  cost)  -  $), 

[(R-k)R-j]  -  (k  -j)  »  (l/|x-X  |J)[Z(t0;  t)  -  Or  sind,  (4.40) 

[(R*k)R*k]  -  (k-k)  =  (l/|x-X|2)[Z(t0;t)  -  C]2  -  1  ,  (4.41) 

CR  x  k)i  =  -(l/|x-X|)(X*j)  *  — (l/|x— X |)r  sinO,  (4.42) 
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(k  X  k)-j  -  (l/jx-X|)((X*l)  -  §]  -  (l/|x-X|Xr  cosd  -  Z)  , 


(4.43) 


(R  x  k)*k  -  0,  (4.44) 

where 

|x-X|l  -  |x|2  ♦  |X|*  -  2x-X  -  |x|2  -  r2  ♦  Z2(t0;t)  -  2[fy  cosd  +  £Z(t0;t)] 

-  |x|2  -  C1  ♦  r2  ♦  CZ(t0; t)  -  CP  -  2^r  cost 

-  ♦  r2  -  7fy  co sd  *  [Z(t0;t)-Cf. 

(4.45) 

Denoting 

Z.(t0;t)  -  Z(to;  t)  -  C  (4.46) 

so  that 


W(t0;t;)  -  -3(l/c)V(t0;t/)Z.(t0;t,y|x-X|  (4.47) 

we  then  have 
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a  x 


E.(x,  t)-i  «  K  J dtcG  '(Ggocng  J dr  r  JddxD(to;  1 0Z.(to;  t  +W(tc;  t  ')](r  cosd-^)/ |x-X  | J,  (4.48) 

0  0  -* 

t.  as 

E.(x,t)-j  -  K  J dt0G /(Qto)G(iJt0)  Jdr  r  J d0xD(to; t ')Z.(tD; 1 ')[ 1  +W(t0; t ')] r  sind/|x-X|\  (4.49) 


0  -x 


T.  «  * 

E.(x,t)-k  «  K  JdtoG/(Qto)G(f2t0)  Jdr  r  JdfixD(t0;  1 '){ [Z.J(t0;  t  'y  |x-X  |3]  -[l/|x-X  |] )  H  +W(t0;  t ')] , 


o  -* 


(4.50) 


(4.51) 


B„(x,t)*l  -  (K/c)  J dteG /(flto)G(Gt0)  Jdr  r  JddxD(to;  t')[l+W(to;t ')]r  sinfi/|x-X|2, 

0  0  HI 

T.  a  s 

B.(X, t)*j  *  -<K/c)  JdtcG  '(PgGtfig  Jdr  r  Jdd*D(t0;  t 'Ml  +W(g 1 7))(r  cosd-^y |x-X|2,  (4.52) 


0  -* 


and 


B.(x,t)*k  *  0.  (4.53) 

We  next  want  to  perform  the  r-0  integrations,  but  we  are  hindered  by  the  fact  that  t'  is  a 
function  of  r  and  d.  Indeed,  we  have  from  Eqs.  (4.30),  (4.24),  and  (328)  that  t'  satisfies,  for 
t,  e  [0, 1.1, 

if  t'e  Kg 
if  t'e  I-(to) 

(4.54) 

(the  non-negative  square  root  is  assumed).  Defining 

A,(tc; t ')  ■  Z,(t0;t/)/^,  (4.55) 


(l/c){52  +  r2-2£r  cosfi  +[vo(t '  - 1„)  -  (v</7)G(ng(t '  -  t^-Q2}02, 
(l/c)fe2  +r2-22jr  cosfi  +  C2}"2, 
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we  may  rewrite  Eq.  (4.54)  as 


(£/cXl ♦A£)UI{  1  -[2/(1  +A^))(r/!j)cosd  +  [1/(1  ♦A£)](r/5)J}Ifl,  if  t7e  Kg 

- 

( |x |/c)[l -2(y\x |  jHr/€)cosd  +  (S/|x|)W]“.  if  t7e  r<g  . 


If  we  now  restrict  our  attention  only  to  situations  for  which 


al%  <  1  (4.57) 

and  further  demand  only  zeroth  order  accuracy  in  (r/2j)  then  Eq.  (4.54)  becomes 
,  ($/c)[l  ♦  A£(t0;t7)]“,  ift'e  I(g 

t  -  t7  «  '  (4.58) 

|x|/c,  if  t7e  I-(g 

with  solution  l'  *  t' (t0;  x.  t)  (that  solutions  to  Eq.  (4.58)  exist  and  are  unique  for  fixed 
tD  €  [O.TJ,  x  «  >  0,  0,  £  £  0),  and  t  k  0  will  be  demonstrated  in  Appendix  Q  which  is  an 

approximation  to  the  exact  solution  t7  -  t'(t„,r,  fi;  x,  t)  of  Eq.  (434)  but  which  is  independent  of  r 
and  6. 

The  condition  of  Eq.  (4.57)  guarantees  that  |5t|(t0,  r,  d;  x,  t)  ■  |t  7(t0,  r,  ,fr;  x,  t)  -  t  7(t0;  x,  t)  | 
is  small,  but  it  does  not  guarantee  that  this  difference  is  small  enough,  where  "small  enough"  means, 
for  example. 


sup  sup  (vo|fit|(t0,r,«;x,tyZM(t0)]  <  1,  (4.59) 

«,  «  [0.TJ  r  c  t0.«],  «  «  [-«.*) 

this  being  a  sufficient  condition  that  every  sheet  moves  only  a  small  fraction  of  its  entire  trajectory 
during  the  time  It7-  t'|.  When  Eq.  (439)  holds,  then  |Z,(tc;t7)  -  Z,(t0;ll)\/ZM(t0)  and 
I V(to;  1 0  -  V(t0;t/)|/v0  are  very  small  and  IXoCgt7)  -  Xo^  *  )  I  ■  1  for  only  a  very  small 
fraction  of  t0  €  [0,  TJ,  hence  the  integrals  in  Eqs.  (4.48M4.52)  remain  quite  accurate  when  t7  is 
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replaced  with  V .  That  some  such  additional  condition  is  needed  is  clear  since  without  it  the  following 
unreasonable  argument  would  be  valid:  Since  we  are  ultimately  interested  in  E,  ->  ■*■  f  or  the  radiation 
limit  then,  given  any  a  >  0,  consider  only  2;  so  large  that  al%<  1;  hence  any  a  >  0  can  be 
accommodated  by  the  condition  o/£  <  1  !  The  condition  of  Eq.  (4.59)  further  restricts  a  so  that  the 
above  unreasonable  argument  is  in  fact  not  valid.  Indeed,  noting  that  using  t '  is  equivalent  to 
assuming  that  all  charge  elements  dS  on  sheet  t0  have  one  and  the  same  associated  retarded  time, 
namely  that  associated  with  the  charge  element  at  r  =  0,  we  see  that  we  may  approximate 
|8t|(tc,  r,  0;  x,  t)  <:  a/c;  further,  from  Eq.  (322)  we  see  that  Z^tJ  £  1/4  vaT  for  all 
tc  e  [0,  TJ.  Hence  a  sufficient  condition  for  Eq.  (4.59)  to  hold  is  that  v0(a/cy(l/4)vor  <  1,  i.e.. 


a!cT<  1, 


(4.60) 


and  we  require  in  the  sequel  that  a  satisfy  this  condition  as  well  as  that  of  Eq.  (4.57). 

Taken  together,  Eqs.  (4.57)  and  (4.60)  simply  mean  that  the  charge  cloud  has  negligible  radial 
extent  for  retardation  purposes.  (In  the  language  of  optics,  we  are  discarding  phase  differences,  which 
lead  to  interference  effects,  along  the  cross-section  of  the  spot.)  For  this  reason  we  refer  to  t/  as  the 
small  spot  siz/e  retarded  time.  Note  that  we  have  made  no  assumption  about  the  magnitude  of  A. 
(which  is  related  to  the  axial  extern  of  the  charge  cloud). 

We  may  now  move  Xd>  Z*  and  the  numerator  of  W  in  each  of  Eqs.  (4.48)  -  (4.52)  outside  the 
r-6  integrals,  replacing  t7  by  V .  When  we  do  so,  and  further  note  feat  fee  integrands  remaining  under 
the  r-f)  integrals  in  the  expressions  for  E,  •  j  and  B,  •  i  are  odd  functions  of  d,  we  get 
E.(x,t)-j  -  0  «  B.(x,t)-i, 

T,  as 

E.(x,t)-i  -  2K|dtoG,(flt^(«t0^cD(to;l,){Z.(to;t,)Jdrr  Jdd(r  cos$^)/|x-X|3 

0  .  ,  0  0  (4.61] 

-  (3/c)V(t0; t ')Z.2(t0; t ') Jdr  r  fdd(r  cosd-$/|x-X|4}, 

0  0 
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«,  «  i  mu 

B.(x,t)-k  -  2K  f dtcG /(flto)G(Qt0)XD(to; ? ,){Z*(to; T 7)  fdr  r  J(hfKl/|x-X|J)-Jdrr  Jdfi(l/|x-X|) 

0  0  0  0  0 

«  «  a  % 

-  (3/c)V(t0; t ')[Z.3(t0; t 7) Jdr  r  JdfiO/|x-X|*)-Z.(tc:t') Jdr  r  Jdd(l/|x-X|2))}  , 


(4.62) 


and 


B.(x,t)-j  -  -(2K/c)  fdtoG/(Qto)G(ft0)xo(to;l/){  Jdr  r  Jdfi(r  costK)/|x-X|J 

o  0  0  (4.63) 

-  (3/c)V(t0;t/)Z.(t0;f,)|dr  rjdfi(r  cosfi-5y|x-X|3}. 

We  see  that  E.(x,  t)  lies  in  the  plane  determined  by  x  and  k,  while  B.(x,  t)  is  perpendicular  to 
that  plane;  the  polarization  of  the  electromagnetic  field  is  thus  specified. 

We  must  next  address  die  remaining  r-fi  integrations.  While  some  of  these  double  integrals 
can  be  evaluated  analytically  in  closed  form,  the  integrated  results  seem  to  us  to  be  too  complex  for 
our  further  use  here  (i.e.,  integration  over  tD)  and  we  will  not  display  diem.  Rather,  we  once  again 
invoke  die  requirement  o/£  -c  1  of  Eq.  (4.57),  this  time  to  expand  the  denominators  in  Eqs.  (4.61)  - 
(4.63).  We  first  write,  as  in  the  RHS  of  Eq.  (4.56), 

|x-X|  «  5(1+Abin{l  -  [2/(1  +A^)] (r/£)cosfi  +  [  1/(1  +Ab](iW “  (4.64) 


so  that,  for  n  =  1, 2,  3, 4,  and  aJ%  <  1 , 


ix-xr 


-  ^l+A^l  +  [n/(l+A*)]cosfi  (r/^)  -  [n/2(l  +/&]  { 1 4(n+2)/(l  +Ai)]cos2fl}(r/&2 


♦  0(( r/^)3)). 


(4.65) 


Using  these  expressions  in  Eqs.  (4.61)  -  (4.63)  and  performing  a  great  deal  of  algebra,  we  find  the 
double  integral  expressions  in  braces  in  those  three  equations  integrate  to  be,  respectively, 
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-(*f2)STO,{A.(lWty*lfl[l  -(IX^-lXl 
-3(V/c)Ai(l  ^Ai)-J[l-(4X2Ai-lXl  ♦AfraTO,]> . 

-(*/2)5TO2{(l  4^)*«[  1  .(4X10AJ  -4Ai+lXl  ♦Ab*aTOJ] 
-3(V/c)A.(l  ♦AfrHl  HjX3Ai-2Aj«lXl  -AirVS)2]} , 


(4.66) 


(4.67) 


(4.68) 

-3(V/c)A.(UAi)‘*,[l-<4X4^-lXl  ♦A^o/^)2]}, 

coned  to  older  (a/tf  inside  the  braces  of  these  last  three  equations.  In  what  follows  we  will  take 
only  die  leading  term  inside  each  of  toe  brackets  (i.e.,  1)  which,  as  can  be  seen,  gives  results  correct  to 
order  \  these  third  order  results  are  more  than  sufficient  for  our  purpose.  The  fields  then  become 


and 


•• 

E.(x,t)  -  -(w2!^  J dtoG  '(QgCKOgXbft.;  t ') 
o 

x  {[l*/&t0;t,)r3a  -  3[V(t0; t VcjA.fv. T ')£  1  WV^(tw; t /))-» ) [iA.(t0; f ')  +  k] 

T. 

B.(x,t)  -  jOu 2K/c4) / dt0G '<QgG(Og XA; t '){ [1  +A&0; t V 


(4.69) 


(4.70) 


-  3[V(t0;tVc]A.(t0;t,)[l^(t0;t,)]-3ffl) 


where 


A.(tc;t')  -  • 


[v„<t'-g  -  (Vy7xxntoxt'-tom  -  TO.  if  t'e  Kg 


-TO. 


if  t  €  rtg 


(4.71) 


and 
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t'  «f'(t0;x,t). 


(4.72) 


To  proceed  further,  we  need  some  additional  properties  of  the  retarded  time  t;.  Because  we 
win  need,  in  addition  to  these  properties,  several  other  results  concerning  the  retarded  time  V  in  order 
to  demonstrate  the  existence  of  the  asymptotic  radiation  limit,  we  systematically  treat  V  in  a 
mathematically  careful  manner  in  Appendix  C.  In  order  to  get  on  with  the  current  development, 
however,  we  simply  note  that  as  a  consequence  of  Theorem  C.3  we  may  rewrite  die  fields  E,  and  B. 
of  Eqs.  (4.69)  and  (4.70)  as 


Tjd.0 


E.(JM)  -  KKV5)  J  dt„G '(flt^ngui ♦A’o,,; t '(t.; x, t))]’3® 

0 

-  3[V(to;t/(to;x,t)yc]  A.(y,  l  '(to;  x,  t))[l  -A&o; !  '(to;  x,  t))]  }  [i  A.(to;  l  '(y  x,  t»  ♦  k] 

(4.73) 


and 


tao 


B.(x.t)  «  j(K'/c£)  f  dt0G,(£2t0)G(Ot0){[lW^t0;t,(t0;x.t))J-> 

0 

-3[V(t0;  l  '(te;  X,  t)yc]  A,(t„;  I  '(t0;  x,  t))[  1  ♦  A&„;  l  'ft;  x,  t))]'«) 


(4.74) 


where  Tc(x,  t)  is  given  by 


To(x,  t)  ■  ( 


0,  if  <  t-  fx|/c  SI  0 

t-|x|/c,  if  0  <  t  -  Jx|/c  £  T0 

T„,  if  T0  <  t  -  |xj/c  <  tjg 

t,(t-|x|/c),  if  tjj')  £  t  -  |x|/c  <  «> 


(4.75) 


wifi) 
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t,(t-|x|/c)  ■  inf  tj[t-|x|/c]. 


(4.76) 


for 


t^ft-lxj/c]  ■  {t0  e  (O.yit.jg  -  t-Jx|/c}; 


and 


(4.77) 


K'  -  xa2K  -  *(mjec *)(iUnyl(a\JT?  -  x[nN(ny2x]2(e  'ltomtc?(mJeWn.)a \AYf . 

(4.78) 

Figure  5  illustrates  foe  determination  of  T0(x,  t),  for  foe  various  ranges  of  t  -  |x|/c  on  foe  RHS  of 
Eq.  (4.75),  as  implied  by  Theorem  C.3,  while  Hg.  6  illustrates  a  generic  T„(x,  t)  as  a  function  of  t  e 
(—»,  oo)  for  a  fixed  x.  Appendix  C  contains  foe  background  required  to  more  fully  appreciate  these 
figures.  Note  that  if  t,  (t  -  |x|/c)  <  t*(t  —  |x|/c)  (in  case  t  -  |x|/c  £  t^TJ),  where 

t,(t-|x|/c)  -  sup  t*‘lt-|xl/c),  (4.79) 

then  we  must  add  ft,  tj  to  foe  ranges  of  integration  in  Eqs.  (4.73)  and  (4.74)  (fra  elaboration  of  this 
last  point,  see  both  foe  paragraph  following  Eq.  (C.30)  and  foe  beginning  of  the  next  section).  Also 
note  that,  as  expected  (see  Eq.  (4.35)),  we  have  from  Eq.  (4.75)  that 

E.(x,t)  -  0  -  B.(x,  t)  if  t  £  |x  |/c .  (4.35) 

The  final  impediment  to  foe  evaluation  of  these  integrals  appears  to  be  this:  we  do  not  yet 
have  an  explicit  expression  for  t '  (te;  x,  t)  when  to  >  0  (t '  (O,  x,  t)  is  given  explicitly  by  Eq.  (C.3)). 
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Figure  5.  Various  possibilities  for  t  -  |x|/c  (where  (t  -  jx |/c)4 ,  i  =  indicates  five  specific 
choices  for  the  value  of  variable  t  -  |x|/c  );  and  foe  ranges  of  values  of  to  for  which 
(tD,  f  y(to;  (x,  t),))  e  D°[r,  G]  according  to  Theorem  C.3.  The  time  variable  is  denoted  here 
by  x  to  distinguish  it  from  t 
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fj(x,  t)  (x  fixed) 


Figure  6.  The  function  T„(x,  t),  t  e  (—>, «),  for  x  fixed.  Note  that  the  function  plotted  here  is  just  file 
"inverse"  of  die  (two-valued)  relation  plotted  in  Fig.  S.  where  file  upper  curve  of  that 
relation  is  first  "fixed  up”  by  removing  that  part  of  its  domain  where  It  is  not  one-to-one. 
The  chape  of  the  portion  of  the  curve  determined  by  tj  can  be  inferred  from  Fig.  5  or,  more 
formally,  from  the  result  (dt/dtXtMdCuVdtKt)  -  l/(dtJdtxC(t))  for 
t  c  ft»(T0)  +  |xj/c,  »)\{(t}.  Finally,  if  there  are  several,  say  NA  >  1,  intervals  of 
constancy  of  then  there  are  NA  places,  tv  i  =  1„..NA,  where  the  map  t  *♦  T0(x,  t)  has  a 
jump  discontinuity  (if  NA  =  0  then  there  are  no  such  places). 
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Such  an  expression  is  given  by  the  unique  solution  of  Eq.  (4.58);  however,  as  discussed  in  Appendix 
C,  we  choose  not  obtain  this  solution  explicitly  since  it  is  not  needed  to  proceed  to  the  radiation  limit; 
indeed,  die  expressions  for  E,  and  B,  of  Eqs.  (4.73)  and  (4.74)  suffice  for  that  purpose.  Equations 
(4.73)  and  (4.74)  thus  give,  adequately  for  our  purposes,  the  final  representation  of  the  finitely-remote 
electromagnetic  fields  in  terms  of  only  the  given  input  parameters  of  the  pulse  and  surface  and  the 
spacetime  point  (x,  t)  of  interest  The  only  approximations  we  have  made  are  (1)  small  spot  size: 
a/%.  a/cT  <  1 ;  and  (2)  yjc  <  1,  which  allowed  us  to  neglect  terms  of  order  (yjc?  in  the  radiated 
fields.  Expressions  fix  foe  customary  finitely-remote  radiation  quantities  may  be  formed  at  this  point 
from  these  Adds  according  to  Eqs.  (4.16H4.20),  but  they  suffer  from  foe  same  deficiency  as  foe  field 
expressions  themselves,  namely,  foe  lack  of  an  explicit  expression  fix  the  small-spot  retarded  time 
f'(t0;  x,  t).  For  this  reason,  we  defer  treatment  of  the  radiation  quantities  until  alter  we  obtain,  in  foe 
next  section,  the  radiation  limit  for  the  fields. 
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V.  ASYMPTOTIC  FIELDS 


The  finitely-remote  Adds,  given  at  spacetime  point  (x,  t)  by  Eqs.  (4.73)  and  (4.74)  in  case  t  > 
|xj/c  and,  more  explicitly,  by  Eq.  (4.35)  in  case  t  $  |x|/c,  depend  upon  x  =  |xj&  in  two  different 
ways;  hence,  so  do  the  finitely-remote  radiation  quantities  calculated  ft  m  them  via  Eqs.(4.17M4.19). 
The  first  dependence  on  x  is  soldy  through  |x|  and  is  embodied  in  the  circumstance  indicated  above 
dut  the  expressions  foe  the  Adds  depend  upon  the  sign  of  t  -  |x|/c;  this  dependence  is  a  causal  one 
and  cannot  be  dispensed  with.  The  second  dependence  is  manifested  by  first  rewriting  Eqs.  (4.73)  and 
(4.74)  respectively  as 

|x|E.(x, t)  -  -K'(l/sii*)[l/UJa(x, t)  ♦  k/aM(x, t)}  (5.1) 


and 


fx|B.(x,t)  -  (K//c)(l/sii*yai(x,t)j. 


(5.2) 


where  4  €  (0,  */2]  is  the  angle  between  x  «  (£> 0,  £)  and  k  (see  Rg.  4)  and 


A5(t0;t'(t0;x,t)){[l  +A£(t0;t'(to;x.t))n 

0 

-  3[Y(to;l,(to;x,t)yc]A.(t0;!,(to;xft)Ml  +  ^(t0;t/(t0;x.t))]-(,»*W)}; 


(5.3) 


and  by  secondly  noting  that  in  general  /^(x,  t)  indeed  depends  upon  x  through  both  |x|  and  «.  This 
situation  is  unsatisfactory  since  die  angular  density  of  a  (conceptually  idealized)  radiation  field  ought 
to  depend  in  fids  second  sense  only  upon  direction  &  and  not  additionally  upon  |x|;  fids  second 
dependence  upon  x,  insofar  as  it  involves  |x|,  should  be  eliminated.  We  accomplish  this  by  defining. 
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for 


ft  -  (slot,  0,  cost),  t  «  (0,  */2] 


(5.4) 


and 


t  €  (-*■>, »), 


(5.5) 


the  functions 


n!(i,t)  -  Um  [|x|E(|x|ft,  |x|/c  ♦  x)J 


(5.6) 


and 


ii(t,x)  -  lim  [|x|B(|x|ft,  |x|/c  +t)], 
Nl— 


(5.7) 


provided  the  limits  exist  If  they  do,  then  we  define  the  asymptotic  instantaneous  radiated  power 
angular  density  at  x,  |x|  »  0,  in  the  direction  ft,  as 


(8P/84'  )_(x,  t)  «  e0cJftnf(ft,t  -  |x|/c)  x  n!(ft,t  -  lx|/c), 


(5.8) 


as  suggested  by  Eq.  (4.17);  other  asymptotic  radiation  quantities  easily  follow.  Of  course  when  x  £  0 
then,  by  Eq.  (4.35),  the  limits  in  Eqs.  (3.6)  and  (5.7)  are  trivial: 


*!(ft,x)  -0  -x!(ft,x)  if  x  £  0; 


(5.9) 
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hence  we  restrict  ourselves  in  the  sequel  to 


x  >  0 


(5.10) 


unless  otherwise  indicated. 

Before  presenting  the  main  result,  we  address  foe  issue,  alluded  to  briefly  after  Eq.  (4.78),  of 
those  special  pulses  -  examples  of  which  were  given  in  Appendix  A  -  for  which  the  integration  range 
in  Eq.  (5.3)  must  be  extended,  for  special  values  of  x  and  t,  beyond  T0(x,  t)  to  t^t  -  |x|/c).  For  such 
pulses  and  such  (x,  t)  we  have,  for  all  t0  e  fF0(x,  t),  t^t  -  |x|/c)],  that  Z(t0;  t  -  |x|/c)  =  0  so  that 
l /(t0;  x,  t)  =  t  -  |x|/c;  hence  Z(tc;  t'(tc;  x,  t))  =  0  as  well  so  that  A,(tc;  f (tD; x,  t))  ■  =  -con> 

and  1  ♦  A^t„;  t(t0;  x,t))  »  l/sirfy.  Since  also  V(tc;  t  -  |x|/c)  =  -vo  then  V(tD;  t'(t0;  x,  t))  =  -vo 
as  well,  hence  foe  additional  integral  that  must  be  appended  for  these  special  cases  is  simply 

VG-I*W) 

/;,(x,t)«cot^[sin>  ♦3(v(/c)cot<>sin1‘*,4]  J  dtcG '(OgGCQg.  (5.11) 

TjCM) 


Denoting 

x,  if  0  £  x  S  T0 

T„“(x)*  T0,  ifT0<x<trtB(f0) 

if  £  x  <  ~, 


(5.12) 


so  that  T„(x,  t)  =  T0"(t  -  jx|/c)  whenever  t  -  |x|/c  £  0,  we  then  have,  using  Eq.  (4.75), 

v« 

£,(|x|*,  |x|/c  +x)oc  J dt0G/(Qt0)G(nt0)  -  (WOftQHOtJti)  -  G*(OT“(x))]  (5.13) 

frw 

so  /p],(jx|8,  |x|/c  +  x)  is  independent  of  |x| ;  hence 
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lim  £,<1*1*.  |*|/c  ♦  t)  *  (l/2)fl'lcot^[sln2>  ♦  3(v  /c)cot*  sin*  ‘  >] 

x  [G\Ots(x))  -  G :(QT"(t))] 


(5.14) 


and  the  existence  of  the  limit  is  demonstrated  for  this  additional  piece.  For  this  reason,  we  will 
dispense  with  further  consideration  of  these  special  pulses;  i.e.,  in  the  sequel  we  will  only  consider 

t,  €  [0.  T„(x,  t)l.  (5.15) 

For  such  t»  we  have  from  Theorem  C.3  that,  whenever  x  »  t  -  |x|/c  £  0,  then 

f(t0;  x,  t)  e  [t0,t0  T/GiOlJ] 

and 

t  -  |x |/c  e  [t0,t0  ♦  T/G(atJ]. 

Also,  for  future  use,  we  plot  in  Fig.  7  the  function  T“(x),  x  2:  0,  corresponding  to  the  function 
T,(x,  t),  t  £  |x|/c  of  Fig.  6  except  that,  in  consonance  with  the  previous  discussion  in  this  paragraph, 
we  take  NA  «  0  in  the  latter  figure. 

The  next  result  establishes  the  existence  of  the  limits  in  Eqs.  (5.6)  and  (5.7)  when  Equations 
(5.4)  and  (5.10)  govern  $  and  x.  The  proof  of  this  theorem  is  given  in  Appendix  D. 


(5.16) 

(5.17) 
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t 


Figure  7. 


“I - 

U(tJ 


The  function  T"(t),  t  £  0,  corresponding  to  the  function  T„(x,  t),  t  £  |x|/c,  of  Fig.6, 
except  that  here  N4  =  0. 
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Theorem  5.1:  Letts  (sin^,  0,  cos $)  with  0  <$£  n/2  and  x  >  0.  Then 


T"(t) 

lim  [  |x |E,(  |x |i,  |x  J/c  ♦  t)]  «  (i  -  £  taafr)K'sini»  cos*  f  dt0G  '<agG(ag 
t*M-  $ 

x  {1  OfV(t0;fl(t0;*,T)ycJcos*}  (5.18) 


trw 


lim [|x|B,((x|i,  |x|/c  ♦  t)]  -  j(K//c)sin*  dt0G '(Qt^Og 


x  {l  ♦  3[V(t0;t(t0;4,t)ycicos«}, 


(5.19) 


where,  if  to  *  0,  then 


v(t0;tl(t„;*,t)yc 


f  (V</C)[i  -  2r-1G<at0)(T  -  g]. 


if  4  *  ic/2 


(i/cos^)(i  -  qojoo  ♦  (2C2(^yc1lwjG(fttoXT  -  g}“),  if  0  <  ♦  <  * 

(5.20) 


with 


C,W  -  1  -  (v,/c)cos* 


(5.21) 


and 


Cj(4)  -  2r-1(v,/c)cos*. 


(5.22) 
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It  is  shown  in  Appendix  D  that  Eq.  (5.20)  may  be  approximated,  to  first  order  in  v^c,  by 

v(t0;tl(t„;<M)yc  -  (vyc)[i  -  2r1G(nt0xt  -gi,  o  <  *  s  *n.  (5.23) 

We  may  dim  use  this  result  to  write  our  final  expressions  for  the  asymptotic  values  of  ref  and 
namely, 

tc!(*.  t)  *  K'siM»K_(*,TXj  x  *)  (5.24) 

and 

*!(*.  t)  «  (K//c)sin4  *.(*,  t)J,  (5.25) 

where 


KJk,x)  «  7.(1, 0;t)  +  3(vo/c)cos0  [7.(1, 0;t)  -  27.(2,  1;t)] 


(5.26) 


for 


T- 

7.(1,  0;t)  -  Jdt0G/(Qt0)G(nt0)  -  (2n)-1(n/G)G  2(Gt0"(x)) 
0 


(5.27) 


and 

t;w 

7.C,  i;x)  -  r»  J  dt0G '(ogG 2(Qt0xx  -g 

0 

-  (3n)-l(w/0){rltT  -  f;(t)]G3(£lf;(T))  «•  G.  3(T0-(x))) 


(5.28) 
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I 


s 

<?_>)  ■  r-,|dt0G'XQt0), 


I,  M  €  N0 ,  s  £  0; 


(5.29) 


here  t  >  0,  0<  *  £  x/2.  and  */I2  is  die  pulse  width. 
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VI.  RADIATION  QUANTITIES 

Having  obtained  Eqs.  (5.24)  and  (5.25)  as  characterizing  the  asymptotic  electric  and  magnetic 
fields,  we  are  now  prepared  to  obtain  the  (asymptotic)  radiation  quantities  delineated  in  Eqs.  (4.16)  - 
(4.19),  as  well  as  some  additional  ones.  We  reiterate  that,  as  per  Eq.  (5.4),  0  <  4  £  it/2. 

A.  Poynting  Vector 

Guided  by  Eq.  (4.16),  we  define  the  asymptotic  Poynting  vector  by 

S^.(x,  t)  *  e0ca|x|**icf(i.  t  -  |x|/c)  x  i£(&,  t  -  |x|/c).  (6.1) 

Bom  Eqs.  (5.9),  (5.24),  and  (5.25)  we  then  have,  to  first  order  in  vjc. 


e0cK/2 |x | -J  sin2^  t  -  |x  |/c)i, 

0, 


if  t  >  |x|/c 
if  t  £  |x  |/c 


(6.2) 


where  JC’i'.&T)  is  in  fact  Kl(k,x)  minus  its  term  containing  (v,/c)2;  namely 


x)  •  /i(l, 0;  x)  *  6(vt/c)cos0  7.(1, 0;  x)[7.(l, 0;  x)  -  2  7.(2, 1;  x)J,  (6.3) 


and  neglect  of  the  (vjc)2  term  in  Kl  is  justified  as  follows.  From  Eqs.  (5.28)  and  (5.17)  we  have 


TTtt) 


0  <  7.(2,  l;x)  £  r-1  J  dt0G,(£2t0)G2(ftt0)[r/G(nt0)]  *  7.(1, 0;x); 


(6.4) 


since  x  >  0  then,  by  Eq.  (5.12),  T “(x)  >  0  hence  7.(1, 0;  x)  >  0  and 
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y_(2.i;tyy_(i,o;x)  *  1. 


(6.5) 


Therefore 


9cosa#v(/c)Ji/.(1.0;t)  -  2  7.(2, 1;  t)|s/72(1.  0;  t)  <  81  iyjcf 


(6.6) 


so  that  the  term  in  Kl  which  is  second  order  in  vjc  may  be  neglected  relative  to  72(1, 0;  t). 


B.  Radiated  Power 

1.  Angular  Density 

From  Eqs.  (5.8),  (5.9),  (5.24),  and  (5.25)  or,  equivalently,  from  (see  Eq  (4.17)) 


(SP/SHOJx,  t)  -  |xP*  S^Jx,  t). 


(6.7) 


we  find  the  asymptotic  instantaneous  radiated  power  per  unit  solid  angle  (*F)  in  die  direction 


i  *  I  sin  4  ♦  kcos4  *  k. 


(6.8) 


for  |x|  >0, lobe 


(«P/W).(x,t) 


e0cK'Jsin24^.(«,t  -  |x|/c), 

0, 


if  t  >  |x|/c 
if  t  £  jx|/c 


(6.9) 


(to  first  order  in  vjc).  We  also  find  using  Rg.  5  and  Eqs.  (6.3)  and  ( 527)  that  the  peak  (in  time) 
radiated  power  per  steradian  at  any  +  e  (0,  x/2]  is  given  approximately  (neglecting  die  term  of  order 
v</c)  by 
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(SP/mr*  (t)  -  e0cK/Jsiii24(2K)-J(ic/Q)2G4(flff0) 


(6.10) 


at  time  t«  |x|/c  +  T„. 

2.  Angularly  Integrated  Power 

Noting  that  4  is  in  fact  the  polar  angle  measured  from  the  (positive)  2-axis,  we  compute  the 
integral  of  (SP/5T )_(x,  t)  over  the  forward  hemisphere  ♦  t)1  ♦  £*  «  |x|l  »  0,  C  2  0,  as  (since 
(5P/54*)_  is  only  defined  for  4  >  0 ) 

tn 

P„(|x|,t)  -  lim[2x  (<*'  sin4/(8P/SH< )_(*, t)].  (6.11) 

4-*0*  \ 

We  find 

• 

eocK^i.(t-|x|/c), 

P.(!*|.t)  - 

0. 

(to  first  order  in  \J€)  where 

-  *{(4/3)/2(1,0;t)  ♦  3(vyC)/_(l, 0; t)[7j(l, 0; t)  -  2/_((2,  l;t))}.  (6.13) 

C  Radiated  Energy 

1.  Angular  Density 

Horn  Eqs.  (4.18)  and  (6.9),  the  asymptotic  total  radiated  energy  per  unit  solid  angle  (to  first 
aider  in  vjc)  is  given  by 
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(8W/5H»).(i)  -  t'CYf1  sin**  /dTj^JM) 

0 

-  eccK/2  rin^ui  ♦  6(v/c)cos*]  J dx/i(l. 0;  x) 

0 

-  12(v^c)co^|dt/.(l,0;x)J.ai;t)}. 

Hie  first  integral  in  the  above  equation  may  be  evaluated  by  writing 

.  -  fr»  t 

Jdt/i(1.0;t)  -  jdx  |  J  dtoG#(Qt0)G(Ot0) 


(6.14) 


(6.15) 


and  then  observing  that  the  integration  domain  for  the  above  triple  integral  is  that  illustrated  in  Fig.  8, 
so  that  a  change  of  integration  order  gives  (using  Eq.  (3.24)  for  t^J*)) 


■IV.W-WM 


JWi(l,0;t)  -  JdscG /(Qsa)G(fisll)  Jdt„G /(Oto)G(flt0)  J  dt 

-  JdsoG  '(OsJGiOs^  Jdt„G  /(Qt0XKnt0Xmin{t^(so).  ^(tj}  -  max{so.  to}) 

X. 

-  { Jds0G'(Qs0^ns0)Jdt0G/(Qt#)G(Qt0)[77G(Qs0)] 

0  0 

».  T. 

♦  J* dsoG  /(Qso)G(ftse)  J* dt0G  /(Qt0)G(Qt0)[7yG(Qt0)l } 

-  (l/3)fr27GW0). 


(6.16) 


Similarly, 


/dt  7_(l,0;x>/.(2, 1 ; x)  -2.  O^TG3#*,)  ♦  i  G‘27TG(fyo)G_s,(ro)  -  (6.17) 


where  GmJl0)  is  given  by  Eq.  (529)  with  s  =  T0.  Hence  we  have 
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t 


Figure  8.  The  integration  domain  for  die  triple  integrals  of  Eqs.(6.15)  and  (6.16),  corresponding  to  the  f0“(x), 
of  Fig.  7.  Given  (So,  U  e  [0,  TJ2,  the  lower  limit  for  die  t  integration  is  s0  ift0  £  so  or  it  is  ta 
if  £  t„;  Le„  it  is  max{s»  tc).  Likewise,  die  upper  limit  is  min  {^(sj. 
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$WW¥)Jk)  -  ^XcK“7WQ)W*(g  *(**,)  ♦  (v0/c)co8*{GJ(flf0) 
-12[G(Cl^s(rj-G^4(!f^])). 


(6.18) 


1  '’’^Energy 

The  asymptotic  total  radiated  energy  is  just  the  angularly  integrated  asymptotic  (energy) 

angular  density,  namely 
to 

Wm  •  lim[2x  f  d*'  fto4'(8W/&¥  )_(*)] 

♦-*«•  $ 

-  ^UO?$G\CBi}  *  ^yJc){G\OiJ  -  UKKarjG^J  -  G^])). 

(6.19) 


3.  Efficiency 

Noting  that  the  photon  pulse  input  energy  for  photons  of  frequency  v  is  (using  Eq.  (2.7)) 


Wta  -  hvna’A  JffOdt  -  hv(x/Cl)a ’AGN^KKflff,,), 


% 


(6.20) 


which  will  be  less  than  the  total  energy  in  the  photon  pulse  if  To  <  x/fit,  we  may  calculate  an 
asymptotic  radiation  efficiency,  using  K7  from  Eq.  (4.78)  and  T  from  Eq.  (3.25).  as 


eT-  ■  W_/Wta  -  (imtfX'/JcymJie  2fejn/:)^)^i(/QftQN(a)]2a  2(AY)2Y 
x  f*G  W„)  ♦  -  12{G_  3(ro)-[G.  4(ToyG(af(>)]})]. 


(621) 


Note  that  (eAYf  is  the  square  of  the  emission  current  density  and,  from  Eq.  (2.6), 
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(6.22) 


■ 


QN(fl)  -  Jg(s)ds. 


Also,  since  for  \l' t  p  t  0  we  have 


*  f 

\G^<iyG Wjl  s;  r-1  J'<h0[G>vnt0yG‘xG?0)]G^£it0)  <  r*1  Jdt0  -  yr 


(6.23) 


then  (noting  G(flT0)  >  0  since  Tc  >  0  by  Theorem  3.2(1)) 

4|G2(flfJ  -  lUG^Gj  -  [G.i4([oyG(ao)])l  /  1  £  1[1  * 24(t,/7)] 

(6.24) 


so  that  tiie  yjc  term  inside  foe  large  brackets  in  Eq.  (6.21)  may  be  ignored  whenever 

\Jc  •  (3/8)[l  +  24 QiJT)]  <  0.05  (0.05  is  arbitrary  but  about  right),  i.e.,  whenever  T a/T  <  1/2.  A 

sufficient  condition  for  tills  is 

1  i»  2(n/Q)T*1  -  «QN(QKco|/U)2  (6.25) 

and  in  that  case  we  have 

e.Td  -  (1/1  &KtXvJc)mt(e  2/eomec)J(hv)-1(ic/G)2[nN(G)]2fl  \AYf  YG  (6.26) 

D.  Spectral  Intensity 


1.  Angular  Density 

We  may  calculate  the  asymptotic  spectral  intensity  angular  density  (52I/5'i#  5 <n)  Ji,  to), 


specified  unambiguously  by  (see  Eq.  (4.21)) 


(5W/8*F  )_(£) 


[do#  (5*1/84'  5oj)_(4,  o>) , 


ftm n  (see  Eq.  (4.19)) 


OT&F  &»)..(*,  to)  -  2£0c1Re(i  •  i&(£,  co)  x  (*V(«,  to)] 


(6.27) 


(6.28) 


(where  a  denotes  Fourier  transform  and  *  denotes  complex  conjugate).  That  this  specification  of  the 
asymptotic  of  spectral  density  is  indeed  the  correct  one  follows  from 

-  fdr(SP/S¥)Jt,T) 

•  •  Jdt  nf(i,t)  x  *i(£,x) 

-  zocH.  •  Jdtl(2it)-W  Jdta  J&(£,  ©)e  ^  x  (2 *)'“  Jdco'i&tf,  a/)e  ^j  (6.29) 

«  -Jdco  i&Xw)  x  (n^’(i.a>) 

-  Jdto  [2eoc  •*&(£,  co)  x  (jOo*(£,  ©)]. 

0 

To  compute  the  spectral  intensity  via  Eq.  (6.28)  we  first  supplement  the  definition  of  £.(£.  x) 
for  x  >  0  in  Eq.  (5.26)  with 

*■_(£,  x)  -  0  if  x  £  0  (6.30) 

so  that  Eqs.  (5.24)  and  (5.25)  bold  even  when  x  &  0  and  are  in  agreement  with  Eq.  (5.9).  Now 


j &(*,  00)  *  (J  x  *)  K'sin*  KJ&  (D) 


and 


i&(«.co)  -  j(KVc)sind  KJt.a) 

with 

m  • 

tf^(*.a»  -  (2*rw  fdt  e  *•**.(£,  t)  -  (2ic)-w  fdt  e  *'*_(£,  t) 

i  o 

•»  • 

*  (2nym { [  1  +3(vo/c)cos$]  J" dt  e^J.(l,0;t)  -  6(vo/c)cos4»  Jdt  e 

0  0 


Using  i  from  Eq.  (6.8)  we  then  have 


(IPV&V  5oj)_(£,  o»  -  2£0cK'W*  \KJ£,  co)  |J. 


The  integrals  in  may  be  performed  as  earlier  (see  Eqs.  (6.15)  -  (6.17)): 


TTW 


JdTe*y_(l,0;i)  *  Jdxe^  J  dt  G  '(flt^KKttg 
0  0  0 

*.  K.V 

-  fdt0G'(QtJG<ntJ  fd re* 

0  I. 

T. 

-  -ioo*1  J dt 0G  '(flt^ftt^  ‘“•(e  *-7XKQ**)-l], 


(6.31) 


(6.32) 


L(2,l;t)}. 

(6.33) 


(6.34) 


(6.35) 


where  this  integral  is  proper  since  both 


(6.36) 


llm  G(Qiy^0<fkJ  •  0 
•.-♦o' 


and  G'(ftto)  »  gCQt^QNCfl)  hold;  and.  similarly. 


jdte‘-7_(2.1;T)  -  -ico-T'1  Jdt0G  e^e"***-*-!]. 


(6.37) 


Hence  we  finally  have 

(S’l/S'F  5co)_(£,  w)  «  Ji'JeocK/2(n/n)3(ft/a})2sin2<> 

x  | Jc!t0G '(QtJGidtJe ““•[e‘“n0(tt*)-l]{l  -3(v,/c)cos<>[l  -2(t/D))\2. 

0 

(6.38) 


2.  Angularly  Integrated  Spectral  Intensity 
Hie  general  expression  for 

an 

(5l/8o>)Jco)  «  lira  [2n  f  c^'sinWl/S'i'  5o>X*.  co)L  (6.39) 

while  straightforward  to  compute,  is  cumbersome  since  cos  9  occurs  inside  lAT^Cfc,  to)  |2;  we  will  not 
display  it  here. 


i 
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VII.  EXAMPLES 


In  this  final  section,  we  present  illustrations  of  our  general  formalism  for  the  following  pulses: 
constant,  linear  ramp,  triangular,  parabolic,  and  sin2.  While  we  have  complete,  detailed  results  for  all 
these  pulses,  we  present  such  detail  here  only  for  the  constant  pulse;  for  the  others,  we  discuss  merely 
those  of  their  features  which  differ  significantly  from  those  of  the  constant  pulse. 

A.  Constant 

A  constant  (or  "flattop")  pulse  is  given  by 

/(t)  =  A,  t  e  [0.  n/Q]  (7.1) 

for  arbitrary  A  >  0.  We  then  have,  for  s  =  f2/t  e  [0,  it], 

g(s)  -  1,  G(Ot)  -  t/(ic/f2),  G'(Qt)  *=  1/it,  ON(fJ)  =  it;  (7.2) 


and 


T  =  2/(it/Q)c^  =  2k-,(vj/c)/(h/G>4Y 


(7.3) 


where 


k  ■  eVe^c. 
Now  for  tc  €  [0,  it/fi]. 


(7.4) 
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WU  *  -  [77(*/a)]ity(j^Q)]-1}  » tc  ♦  (2 /<#C 


so  ^  has  unique  local  minimum  at 

tA-jB(7)  -  (n/a)[r/(Jt/fl)]1/2  -  V2/<bp  -  ^2k-w(vJc)m/(AY)m. 
We  then  see  that  tft-ta(7)  <  x/Q  iff  T  <  k/Q  so 


!  K*fn- 

if  S 

n/Cl, 

if  W 

where,  observing  that  To  depends  upon  (or  7)  which  in  turn  varies  directly  as  Am 
A),  we  have  distinguished  "strong  pulse”  (S)  and  "weak  pulse”  (W)  cases  as  follows: 


S:  Wj/Q  >  J 2/n  or  T  <  x/ft 
W:  (Oj/Q  $  yJlI-K  or  T  £  it f£l . 


The  situation  is  illustrated  in  Fig.  9.  Continuing,  we  find 


uu 


2*^(7). 

if  S 

n/a  ♦  T, 

if  W 

and,  for  t  *  ^(TJ, 


C(0  -  0/2){t  -  [t2  -  4(n/Q)TJU2}  =  tfi)  -  tjd). 


(7.5) 


(7.6) 


(7.7) 


(or  inversely  as 


(7.8) 


(7.9) 


(7.10) 
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UO  =  t.  +  [V(1C/Q)K1 


Figure  9.  The  strong  (S)  and  weak  (W)  pulse  cases  for  the  flattop  pulse  of  Section  VILA. 
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Then  T"  is  explicitly  specified  for  the  strong  and  weak  cases  by  using  Eqs.  (7.7)  and  (7.10)  in  Eq. 

(5.12);  further 


7.(1, 0;  x)  -  (27r)-1(Jt/n)[t0‘(ty()i/fl)]2,  (7.11) 

7.(2, 1; x)  =  (3n)-‘(it/ft)[t(r(xy(it/fl)]JT-,[t  -  i-T^x)],  (7.12) 


and 


A^(*,x)  -  (2it)*J(iKfl)1rro'(ty(Ji/0)r(l  ♦  6(v,/c)cos${l  -^(x^n/fDJT'^x  -iTjXx)]}) ; 

(7.13) 


also. 


K'  -  (l/4Jm?{mJeMa)a\AYf 


(7.14) 


and 


eocK/2  «  ( l/16)ir2icJin(l(n/0)Ja  4(A Y)4 .  (7.15) 

We  are  now  prepared  to  display  the  radiation  quantities.  If  t  >  |x|/c  andx  =  t  -  |x|/c  then 
the  radiated  power  angular  density  is 

(SP/8'F)_(x,t)  -  ^Khn.sitf^aVY/ft^x)]^  +6(vyc)cos0{l  -|T0"(x)[x  -  if 0"(x)y(x/Q)7l) 

-  J-iem.sin^  fl%4Ymo-(x)r{  1  -  4kcos^(AY)T“(x)[x  -lf0’(x)]  ♦  6(v,/c)cos*} 

(7.16) 
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where  we  have  used  Eq.  (7.3)  to  obtain  the  second  equality  in  Eq.  (7.16).  The  approximate  peak  (in 
time)  power  (gotten  by  neglecting  the  vjc  term  in  the  first  line  of  Eq.  (7.16))  is 

[77(jc/Q)]2.  if  S 

(SP/S'i'rV*)  -  ( l/64)ic3inesin24  a  x  (7.17) 

1,  if  W; 

note  that  in  the  strong  case  this  may  be  rewritten,  using  Eq.  (7.3),  as 

(SP/5'F)r*k(t)  -  (l/16)ian,sin>  (\Jcfa*(AY)2  (S).  (7.18) 

The  angularly  integrated  power  and  approximate  peak  power  are,  respectively, 

P„(|*|.t)  -  (l/64)icKhn„fl\AY)4rf0“(x)l4(i  +  3(v(/c){l 

(7.19) 

and 

’ 

[r/K/fl)]2,  if  s 

pr*(|x|)  -  (l/48)KKhn/i4(AY)4(it«i)4  x  '  (7.20) 

1,  if  W 

where,  in  the  strong  case,  the  last  may  be  rewritten  as 

Pr*(|x|)  -  (l/12)manc(vt/c)2fl4(AY)2  (S).  (7.21) 

The  total  radiated  energy  angular  density  is 
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(8W/ST)J£)  -  -iKhn^n2*  a4(AY)4(ji/n)7'  T*{  1  ♦  (v,/c)co*[l  -jT^n^DTl} 

-  ^riAn,sinJ<Kv</c)fl\AY)Jro3{l  - .1kco8«AY)T’  ♦  (v,/c)coaM 

2VJK1/2(vo/c)V3(AY)V5[l  +  |(v/c)cosfl,  if  S 

« 

x2(v(/c)(AY)i(ji/n)J[l  - ^KCoaKx/Q^AY  ♦  (vo/c)coa>],  if  W 


and  the  total  radiated  energy  is 


W_  -  ^inc2m.(v(/c)fl4(AY)JT3[4  -^tcAYT’+lOtfc)] 


r2  ^  l  j 


^3tmea4  x 


23/ikw(v0/c)«(AY)vj[4  «•  ^vjc)}. 


*?(vJc)(AYY(xKin±  -^k(x/Q)2AY  +\(yjc)). 
To  compute  (he  radiation  efficiency  we  note  that 


if  S  (7.23) 

if  W. 


Wta  «  Jthva2AT0 


(7.24) 


so 


^.^(vycXhv)*^  \AYfY  fB  [1  -  JLkAYT’  +i(v,/c)) 


^m,(hv)V  x 


2K(vd/c)2(AY)Y[4  *^yjc)), 
x*(vo/cXAY)lY(x/ft)J[i  -^k(k/Q)2AY  +1^)], 


if  S 

if  W  . 
(7.25) 


Finally,  the  spectral  intensity  is  given  by 
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(8tyS*F&B)_(*<0)  -  (l/ie^-’K^.sin^  a4(AYf[M(ai)/G?}2 


(7.26) 


where 

A/(©)  •  |c,(©)  -  (coscoT0  ♦  ©T0sinaT0  -  1) 

♦  3(vyc)cos$  jc^a)  -  (cosqjI0  +  ©T0sinaT0  -  1)  -  k(-^X-jj-){C,2((0) 

-  2©T0cos©T0  -  ((©I,)2  -  2]sin©T0}| 

♦  |s,(to)  -  (sin©Ta  -  ©^cosal,,) 

♦  3  (vyc)cos  $  |SjCO)  -  (sintoTo  -  ©^cosaTy  -  *C_|.)(._p{52(a) 

-  2coT0sinojTo  -  [(©f,,)2  -  2](l  -  a>s  ©t„)}Jj 


«V«Sr. 

ca(©>  « (.fijlr1  fdy  yq  cost/T— (y  ♦  1)1  (q  -  l,  2) 

’  ®P  o  <*V  y 


and 


<VJX 

S,(©)  ■  (fiJOLr1  |dy  y*1  sin(V^^(y  ♦  -)]  (q  -  l.  2). 
*  o>P  {  %  y 


Note  that  In  the  strong  case. 


©L  - 


and  (©y^o  -  1 


(S) 


(7.27) 


(7.28) 


(7.29) 


(7.30) 


©T0-it“  and  ((0ji/2*e  -  ^  (W).  (7.31) 

B.  Linear  Ramp 

Hie  linear  ramp 

At)  •  A[t/(x/a)L  te[0,x/a],  A>  0  (7.32) 

yields  results  very  similar  to  those  for  the  constant  pulse,  with  the  strong  pulse  regime  resulting  when 
T  <  (1/2)  (x/Q)  and  the  weak  pulse  resulting  when  T  t  (1/2 Xn/O). 


C  Triangular 

For  the  triangular  pulse 

( t/(x/Q), 


/(t)  -  2A  x 


1  -  [t/(x/n)), 


if  0  £  t  £  (l/2X*/0) 
if  (1/2X*/Q)  £  t  £  iua 


(A  >  0)  we  find 


T  -  4/(x/Q)a>p  -  4ir,(v,/cy(Ji/Q)AY 


(7.33) 


(7.34) 


and 


tJU  -  (a/Q)  x 


Kt,/(x/Q)]  ♦  {mwwtmjwr1, 

t V(*«J)J  ♦  (77(n/f))](l  -2(1  -[tydi/Q)]}2)'1, 


if  0  £  t0  £  (l/2X«/fl) 

if  (1/2 Kx/Q)  £  t0  S  x/Q. 

(7.35) 
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Additionally, 


1  -  [77(K/f2)][ty(K/Q)]-\ 


if  0  £  tc  £  (1/2)  (Jt/ft) 


1  -4[77(X/Q)]{1  -(ty(ii/Q)]}(l  -2{I  -(ty(*/Q)]}Jp,  if  (1/2)  (x/Q)  <J  t„  <J  x/Q 

(7.36) 


and  it  is  easy  to  check  that  (dH^/dt^tJ  >  0,  unless  te  *  (1/2)  (x/fl)  where  this  second  derivative 
does  not  exist;  since 


( AJA^K/a )  -  1  >  0 


(7.37) 


then  t^tj  has  exactly  one  local  minimum  in  (0,  x/Q),  say  at  tD  =  taBjD(7).  Now  since  the  expression 
for  tike  derivative  in  the  first  line  of  Eq.  (7.36)  is  zero  iff  tD  is  equal  to 

t£L (T)  m  (jt/QXTXOt/O)]18,  (7.38) 

and  €  [0,(l/2Xx/Q)]  iff  77(x/fl)  £  1/8,  then  -  is  given  by  whenever 

77(*/»)  £  1/8; 

t*«h(7)  -  Cih(7)  -  (*/Q)[7T^Q)]w  If  77(x/tl)  £  1/8.  (7.39) 

On  die  other  hand,  if  77(n/Q)  >  1/8  then  (dtBB/dt0X(l/2XJt/^))  <  0  so,  using  Eq.  (7.37),  we  conclude 
that 
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e  «i/2 MO),*/a) 


(7.40) 


and  is  gives  by  (be  appropriate  root,  t^ljT) ,  of  die  quartic  equation 
4{770c/Q)]{i  -  [ty(»n)]}  -  (1  -  2{i  -  [tyoc/Q)])2)1. 


(7.41) 


namely. 


t^t)  -  tfjl)  if  77(n/fl)  >  1/8.  (7.42) 

We  will  display  t^(7)  shortly,  but  we  now  wish  to  make  an  important  point:  for  our 
triangular  pulse,  Eq.  (7.40)  tells  us  that  To  <  Jt/fl,  always;  Le.,  in  our  model  with  cutoff  we  may  never 
encompass  the  full  triangular  pulse  but  only  some  proper  initial  portion  of  it  On  foe  other  hand,  we 
see  from  Eq.  (7.41)  (or  from  Eqs.  (7.46)  -  (7.48)  below)  that 

Mm  d,(7)  «  s/fl  (7.43) 

80  we  can  in  principle  get  as  close  as  we  please  to  foe  lull  pulse  by  taking  77(ic/Q) 

I*  2k-‘(v</c)/(k/Q)2AY]  large  enough.  In  analogy  with  foe  constant  pulse,  we  call  foe  case 
T  £  (l/8X*/£2)  foe  "superstrong”  (SS)  case  (because  here  foe  cutoff  occurs  during  foe  rising  portion 
of  foe  pulse)  and  foe  case  T  >  (1/8)  (n/fl)  the  strong  (S)  case  (here  foe  cutoff  occurs  during  foe  falling 
portion  of  trie  pulse);  there  is  no  weak  case  but  we  also  distinguish  foe  case  of  foe  "weak  limit”  (WL), 
where  7Y(n/Q)  -»  ■».  We  then  have,  in  analogy  with  Eqs.  (7.7)  -  (7.9), 
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<l'Un 

if  T  £  (l/6Xn/Q) 

(SS) 

t?Un 

if  T  >  (1/8XH/0) 

(S) 

(7.44) 

u/n, 

if  77(ic/Q)-><» 

(WL) 

and 


wygUn 


<JU« 


+T{1  -2[1  -  t^OXn/Q)!*}'1, 

(n/fl)  ♦  T, 


if  SS 
if  S 
if  WL 


(7.45) 


and  we  may  proceed  by  cases  to  the  radiation  quantities.  Finally,  as  promised,  we  have 


t£L(7y wo)  -  i  -  {2  -r2  +  2[T/(iuaM‘l)ul) 


(7.46) 


where 


R  *  +(4  +  2-w(i?.  ♦  Aj]w 


(7.47) 


for 


A*  -  ({£  ♦  [77(^0)?)  ±  [77(a/Q)]{  *  ♦  [77(^0)]l)“)'fl.  (7.48) 

D.  Parabolic 

The  parabolic  pulse 

At)  -  4A[t/(jtfQ)J{l  -  It/Ot/O)]},  t  e  [0,  rc/Q],  A  >  0,  (7.49) 
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yields  results  very  similar  to  those  for  the  triangular  pulse.  Once  again  no  weak  pulse  is  possible; 
further,  since  the  parabolic  pulse  is  defined  as  (me  piece,  it  is  not  as  natural  here  to  distinguish 
between  a  superstrong  and  strong  pulse  (although  this  of  course  may  be  done)  as  in  the  triangular  case; 
we  may  thus  consider  only  a  strong  pulse  and  the  weak  limit  Also,  we  have 

(dt ^/dgco  «  1  -  6[77(*/0)]{l  -  [tyot/fl)]}  /  [tyOi/O)]^  ♦  2{1  -  IVOt/O)]})2  (7.50) 

with 

T  -  3/(x/fl)o^  (7.51) 

and  (d^/dt’XU  >  0  on  [0,  xJQ]  so  determination  of  unique  involves  solving  the  quintic  in 

to/(x/fl)  derived  from  Eq.  (7.50)  and  this  cannot  be  done  analytically  in  general  -  it  must  be  done 
numerically. 

E.  Sine-squared 
The  pulse 

J( t)  -  AshrMwIt/Ot/Q)]}  =  Asin2(Qt),  t  e  [0,it/Q],  A  >  0,  (7.52) 

yields  results  similar  to  those  for  the  parabolic  pulse,  exhibiting  only  a  strong  pulse  and  a  weak  limit 
To  find  ^^(7)  we  must  solve 

stant,/l2nt0  -  sin2£ltj  *  1/2^2  n[T/(n/a)]U2  (7.53) 
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(x-sinx  =  0atx  =  0  and  (d/dxXx  -  sin  x)  =  1  -  cos  x  >  0  on  (0,  K)  so  x  -  sin  x  >  0  there);  this 


has  unique  solution  since  sin£)tl/[2at0  -  sin  2fltJ  strictly  decreases  from  » to  0  on  [0,  k/£1]  j  to, 
but  of  course  it  must  be  found  numerically  in  general. 
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APPENDIX  A:  A  SPECIAL  CLASS  OF  PULSES 
In  this  appendix  we  exhibit  a  class  of  pulses  each  member  of  which  has  the  property  that  its 
associated  map  t0  *-»  t^tj,  tD  e  [0,  Tj,  possesses  a  positive-length  interval  of  constancy  and, 
furthermore,  ?„  =  n/Q  (the  full  pulse  width).  Classes  with  these  properties  other  than  the  one  we  will 
exhibit  also  exist  but  we  will  present  only  one  such  class.  In  particular,  we  will  prove  the  following. 

Theorem  A.1:  Fix  Q,  v*  Y  >  0. 

(i)  If  T  £  r IQ,  0  <  a,  <  2/3,  and  a,  <  ctj  <  r/OT  then  the  pulse 

/^(s)  *  2flJ(l/IeJ(Y/v0ymcEj}pOiOj(s)  .  s  e  [0, it] .  (A.1) 

where 

tl/(my(l  -a^cgt^Ul  -2a,)/a,QTIs  «•  [2  -3a,]},  if  s  e  [0,0,011 
P<w(s)  *  l/(OT  -  s)2  ,  if  s  e  [ctjOT,  a,OT]  (A.2) 

[l/(OT)3(l  —  Oj)2o,]s  ,  if  s  e  [oijOT, n], 

has  the  property 

-  T  for  ah  t0  e  [OjT,  c^T].  (A.3) 

(ii)  If,  in  addition,  T  ^  3(r/0)  and  a,  <  ji/OT  (<  1/3)  then  there  exists  o£  e  (a,,  x/OT)  such 
that  if  a,  e  (Oj,  r/OT)  then  (a)  Tc  =  TJT,  G^]  =  r/O,  where  T  =  2/c^Naaj(0)  (see  Eqs.  (3.25)  and 
(3.36));  and  (b)  t^tf.)  <  T. 

Proof,  (i)  It  is  easy  to  check  that  p^^  is  continuous  on  [0,  r]  and  also  (strictly)  positive 
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there  (since  o,  <  2/3).  Next,  denoting 


A  -  max  /,  (s) 


(A.4) 


we  rewrite  as 


LA s)  -  A{2fl2/le3(AY/voym.eJ}paiaj(s)  «  Ag^s) 


(A.5) 


and,  noting  that 


max  g0  „(s)  «  1, 


(A.6) 


we  then  have 


g^Cs)  =  2(02/<flfow(s). 


(A.7) 


Noticing  from  Eq.  (3.2S)  that 


2(tf/a^)  «  on^o)  •  nr 


(A.8) 


and  from  Eqs.  (2.4)  and  (2.6)  that 


■ 

-  /  g^U')*'.  S  €  [O.K], 


(A.9) 


we  find  from  Eqs.  (A.1)  and  (A.7)  that 
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[WJTftl-a^ogt-ta  -  2a,)/a,flT]s 1  [2  -3a,]s}, 


G^sVQT  - 


l/(OT-s), 


(l/2)[l/flT(l  -aflllVaJUUIft*  *  [2  -30,]}  , 


if  s  e  (0,  a,  £2T| 
if  s  e  [Oi«T,  otjOT] 
if  s  e  [a,flT,  it] . 

(A.  10) 


Hence,  using 


OtJ.0  »  Ot0  +  TOOtJ], 


(A.  11) 


we  find 


^W'o)  * 


flto  ♦  (flTftl  -  a^a.n  -[(1  -2a,)/oinT](Gt0)i  *12-  30,1(00},  if  ta  e  [0,  a,T] 
OT,  if  t„  e  [a,T,  o,T] 

nt„  ♦  2(OTKl  -Oj^/Ul/OjCnTfKQto)2  +  [2  -  3a,]}  ,  if  t„  e  [aJ.n/Q] 

(A.  12) 


so  Eq.  (A. 3)  holds.  (Note  that  t^^CO)  -  ~,  as  required,  but  that  the  denominators  in  the  first  and 
third  lines  of  Eq.  (A.12)  are  never  0  since  a,,  a,  e  (0, 2/3).) 


(ii)  Let  T  £  3(tt/0).  To  prove  (a),  we  will  show  that  if  0  <  a,  <  jt/OT  then  there 
exists  Oj’  €  (a,,  u/QT)  such  that  if  a,  e  (a,*,  n/OT)  then  (dtmvOiOi/dt0)(t0)  £  0  for  all  t0  e  (0,  x/fl]  so 
that  To  =  jt/O,  by  Eq.  (3.36);  note  that  for  such  a,,  a,  the  conditions  of  (i)  hold  so  that  enjoys 
file  property  given  by  Eq.  (A.3). 

We  begin  by  noting  that,  for  0  <  a,  <  2/3  and  o,  <  a,  <  it/flT,  t^  «,«,(to)  ^  differentiable 
on  (0.  x/Q]  since  g^  is  continuous  there.  The  stationary  points  of  t^^  are  thus  found  from 
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-  i  -  nro^sVG^s)  -  o,  «t0  -  s  e  (o,  it). 


(A.  13) 


where,  from  Eq.  (A.10), 


oro^syG^s) 


(flnyo-o,)2^  -2[(1  -2<xiya1GT]s  ♦  [2  -3a,]} 
*  s2{  -K1  -2a,ya,Qt]s  ♦  [2  -  3a,]}2, 

1, 


41(1  -a^/a^s/m/a^Gri^s2  ♦  [2 -3a,]}2  , 


if  s  e  (0,  a,GT] 
if  s  e  [a,GT,  a,GT] 
if  s  e  [o^GT,it]. 

(A.  14) 


If  s  e  (0,  a,  GT]  then  Eq.  (A.13)  becomes,  using 

o  ■  s/a,GT  (0£©£  1), 


(A.15) 


SRo)  •  (1  -2a,)20‘  -  2(1  -2a,K2  -30,)©*  *  (2  -  3a,)2©2  *  (2/a,Xl  -a,)2a  -2 a,)© 
-  (l/a,Xl  -  0,^(2  -  3a,)  -  0  -  (o  - 1)^0) 


(A.16) 


for 


S^a)  m  (l  -2a,)*0s  -  (1  -2a,X3  -4a,)©3  ♦  (1  -a,)2©  +  (l/a,Xl  -a,)2(2  -3a,). 

(A.  17) 


Hence  ©  *  1,  i.e.,  t„  =  a,T,  gives  a  stationary  point  of  ^^(tj.  Further,  as  we  demonstrate  in 
Lemma  A 2  below,  4g(0)  has  no  roots  in  (0, 1]  when  0  <  a,  <  1/3  so  that  o  =  1  is  the  only  root  of 
JJo)  in  10, 1);  i.e.,  tc  *  a,T  is  the  only  stationary  point  of  tmaA  in  [0,  k/G],  Since  dt^^/dt,,  is 
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then  we  have 


continuous  on  (0,  o.T]  and,  from  Eqs.  (A.13)  and  (A.14),  lin^dt^  a  =  —» 

i  ,-*<r  n 

(<K,%^oXV  <  0  for  tc  €  (0,  a,T)  (0  <a,  <  1/3).  (A.18) 

Next,  if  s  e  [otjCIT,  04  QT]  then  Eq.  (A.13)  becomes  1  =  1,  which  is  satisfied  for  all  s  e  [a,flT, 
OjflT],  i.e.,  for  all  t0  e  la, T,  04 T];  of  course  this  also  follows  directly  from  Eq.  (A.12).  So  we 
have 

(dUo.o/rtoXg  *  0  for  t„  e  I04T,  04T].  (A.19) 

Lastly,  If  s  e  [a,f2T,  k]  then,  using 

d - S/C4OT  (l^d^Ji/otjflT),  (A.20) 

Eq.  (A.13)  becomes 

H«*)  ■  afa4  *  204(2  -  304)^  -  4(1  -  04)^  +  (2  -  304)2  -  0  *  (0  -  1)HW  (A.21) 

far 

!^(d)  ■  ot^d3  ♦  ajd2  ♦  04(4  -  50j)d  -  (2  -  3a,)2.  (A.22) 

Hence  d  =  1,  i.e.,  tc  =  04 T,  gives  a  stationary  point  of  ^^(t,,),  in  agreement  with  Eq.  (A.19). 
Further,  since 

(dJVdd)  (1) »  — 12£(l/3)  -  04)3(1  -  04)  (>L23) 
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then  6*1  gives  a  local  maximum  iff  ctj  e  (0, 1/3).  We  now  rewrite  ^Kj(d),  using 

2  ■  A  - 1  [0*2*  (a/OjQT)  -  U.  (A-24) 

to  get 

(l/c£)!)^  (A)  -  (A  - 1)*  ♦  4(6  - 1)*  ♦  (4/a,)  (A  - 1)  -  4[(l/a,)  -  l][(l/otj)  -  3] 

*  2*  ♦  42*  ♦  (4/<Xj)2  -  4[(l/Oj)  - 1]  [(1/a,)  -  3]  *  ^  (D  (AJ5) 

and  seek  the  zeros  of  9^,.  Using  standard  techniques  for  cubic  equations  we  find  for  Eq.  (A. 25)  die 

ftisffimlnant 

A(a,)  «  (4/3)3[(l/o,)  -  (4/3)]3  ♦  4{[(l/a,)  -  (4/3)]*  ♦  (1/27)}*  >  0  ,  (A26) 

where  the  last  inequality  holds  for  a  <  2/3.  Hence  has  only  one  real  zero  in  (-»,  «*»)  and  we 
denote  die  value  of  2  that  makes  ^,(2)  =  0  by  2,(a ,).  In  fact, 

20(a,)  -  ^[(JL)  -(!)]*  *(^)}  ♦  A^J*  ♦  ^{[(-L)  -(!)]*  +(J-)}  -  A^a,)]*-  (1) 

(AJt7) 

but  we  do  not  need  this  explicit  expression  for  our  proof.  (However,  it  is  of  practical  use  ~  see  later.) 
We  will  demonstrate  in  Lemma  A.3  below  that 

2„(a,)  >  0  for  0  <  a,  £  n/OT;  (A2S) 
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we  then  have 


ctj  ♦  041,(04)  >  Oj  for  0  <  a,  £  it/fiT  (A.29) 

so  Out,  in  particular, 

K/OT  *  WttVLJlit/ttn  >  ji/QT.  (A.30) 

Since  04  ►*  a,  +  042^(04)  is  continuous  on  [0,  x/ffl]  then  there  exists  04  e  (a,,  x/QT)  such  that 

°4  +  >  **01  for  04  e  (04,ic/QD  ,  (A.31) 

so  that 

Io(04)  >  X/04OT  -  1  for  04  e  (04,  x/flT);  (A.32) 

hence 


*.(04)  *  2.(04)  ♦  1  >  n/OjQT  for  04  e  (04,  x/fiT) 


(A.33) 


and  so  d  *  1  is  the  only  root  of  IK,  in  [1,  Jc/OjfJT].  Thus  tD  *  04  T  is  the  only  stationary  point  of 

W«a(0  ta  (®iT*  Port11®*  since  6  =  1  gives  a  local  max  for  H  (when  04  <  1/3)  and  HO)  = 
0,  then 


H(6)  <01brde(l,  X/04OT] 


(A.34) 
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since  is  continuous  on  [1,  x/a,QT].  But 


<<*U*-**.XO  •  H^TKc^t./o,!?  -  (30,-2)]*.  tc  e  [o,T,x/Q]  (A.35) 

so  finally  we  have,  using  Eq.  (A.  34), 

<0  for  tc  e  (o,T,  x/Q]  (o,*<a,<  x/QT).  (A.36) 

Taken  together,  Eqs.  (A.  18),  (A.  19),  and  (A.36)  show  Te  =  x/Q,  so  (a)  is  proved 
To  show  (b),  we  first  note  that  since,  by  file  above,  T0  =x/Q,  and  since 

^^(x/Q)  «  x/Q  *  T  -  x/Q  «■  2/cojN.^Q)  (A.37) 

then  we  must  show  that 

x/Q  +  Z/o^N^Q)  <  T.  (A.38) 

Computing 

K 

N^Q)  -  (1/Q)  Jgj^sjds  -  [l/tojT(l  -  a,)1]  [2  -  3a,  +  (l/c^Xn/QT)2) ,  (A.39) 

0 

Eq.  (A.37)  becomes,  for  a,  <  2/3, 

-c£  ♦  (1/2)11  ♦  3(x/QT)]aj  -  (x/QDo,  ♦  (l^Xx/QTftl  -(x/QT)]  >  0.  (A.40) 
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We  will  demonstrate  below  in  I-emma  A.4  that  this  last  equation  indeed  holds  if  * /QT  <,  1/3  and  a,  < 
x/flT  so  that  (b)  is  proved.  ■ 


We  now  present  the  three  lemmas  referred  to  in  die  proof  of  the  theorem. 

Lemma A.2 :  Let0<a<l/3andforael define 

F«(o)  ■  (1  -2a)Jo3  -  (1  -2aX3  -4 arf  *  (1  -afa  *  (1/aXl -a)*(2  -3a). 

Then  F«(o)  >  0  for  o  e  [0, 1]. 

Proof.  From  (l/3)a(4  -  5a)  >  0  it  follows  that 

(1  -a?  «  1  -2a -a2  >  1  -(10/3)0  ♦  (S/3)a?  -  (8/3)[(3/8)  -  (5/4)a  ♦a2] 

-  (8/3)  [(1/2)  -  a]  [(3/4)  -  a]; 

noting  (2/3  -  a)  >  a  we  then  have 

(2  -3aXl  -  a)2  >  a(l  -2aX3  -4a). 

Hence,  for  a  e  {  0, 1], 

(l/aX2  -3aXl  -  a)2  >  (1  -2aX3  -4a)  £  (I  -2aX3  -4a)o 


(A.41) 


(A.42) 


(A.43) 


(A.44) 


so 


(1  -2a)lo3  ♦  (1  -a)*o  +  (1/aXl  -a)*(2  -3a)  £  (1/aXl  -a)*(2  -3a) 
>  (1  -2a) (3  -4a)o  *  (1  -2a) (3  -4a)o2  , 


(A.45) 


from  which  it  follows  immediately  that  Fa(c)  >  0.  ■ 


Lemma  A.3 :  Let0<a<  1/3  and  for  X  e  1  define 

H.(X)  ■  X3  +  4X2  ♦  (4/o)I  -  4£(l/a)  -  l][(l/a)  -3]. 

Let  X0(a)  denote  the  unique  zero  of  Ha  (as  per  Eq.  (A26)  ff.).  Then  X0(a)  >  0. 

Proof.  It  is  sufficient  to  show  that  Ho(X)  <  0  for  all  X  £  0.  To  this  end,  we  let  £ 

SJe  ■  -0^(1)  «  oft?  -  40^  ♦  4oX  ♦  4(1  -o)(l  -3o),  X  e  B 

and  show 

iyX)  >  0  for  all  X-  2:  0. 

If  IT  *  0  then  Eq.  (A.48)  clearly  holds,  so  suppose  X  >  0.  Now  it  is  clear  that 
aX7  -  4aX  +  4  >  0  for  all  X  e  1 

(considered  as  a  parabola  in  X )  so  that 


(A.46) 

-X  and 

(A.47) 

(A.48) 

(A.49) 


ce*3  +  4oX  >  4a3X7  for  X  >  0; 


(A.50) 


tfF  ♦  4al  ♦  4<l-aXl-3a)  -  4a21 2  >  0  for  I  >  0  (A.51) 

and  dlls  is  precisely  Eq.  (A.  48).  ■ 

Lemma  A.4 :  Let  jt/ft T  <,  1/3  and  for  a  e  E  define 

yo)  «  -a3  ♦  (1/2)[1  ♦  3(ii/OT)]a?  -  (x/OTXx  ♦  (l/2X*/fiT)J[l  -  (it/CTT)].  (AJ2) 

Then  J^/a)  >  0  for  a  e  (0,  * SOT). 

Proof.  This  follows  immediately  from 

W°)  *  (1/2Xk/OT)2[1  -(K/flT)]  >  0,  (AJ3) 

J^x/ftT)  «  0,  (AJ4) 

and 


(dl^daXo)  -  -3[(l/3)  -  a][Jt/ftT)  -  a]  (AJ5) 

since  the  derivative  is  negative  for  a  <  min  { 1/3,  it/ftT).  ■ 

As  a  concrete  example  of  a  pulse  in  the  class  of  the  theorem,  we  choose  T  =  4(n /(l),  04  = 
1/32,  and  04  *  1/8;  that  such  an  04  is  adequate  can  be  verified  by  the  fact  that  it  satisfies  Eq.  (A.31), 
with  Eo(0Cj)  given  by  Eq.  (A2T).  We  leave  ft,  Y,  v0,  hence  A  and  eo?,  nonspecific.  We  then  find 
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(gfllft^CXs/x)  +  61], 


8un.u«(s>  -  2(Q/(op)2(l/xJ)  x 


l/[4-(s/x)]2, 

(8/49XS/X), 


if  s  e  [0,x/8] 
if  s  €  [x/8,  x/2] 
if  s  e  [x/2,x]. 


ON^^O)  -  (34/49x)(Q/«p>*. 


<W«0»  »  (49/17)  x 


(8/3 1)*[  -240(s/x)J  *  61(s/x)J, 

if  s  €  [O.x/8] 

1/14  -  (s/x)]. 

if  s  6  [x/8,x/2] 

(l/49)[4(s/x)J  ♦  13], 

if  s  e  [x/2,x]. 

and 


®*w,un.i/»(s) 


x{(s/x)  ♦  (31/8)2[-240(s/x)2  ♦  ei^x)]'1}, 
4x, 

x{(s/x)  ♦  49[4($/x)2  ♦  13]"‘) , 


if  s  6  [O.x/8] 
if  s  e  [x/8,x/2] 
if  s  €  [x/2,  x]. 


(A.56) 


(AJ7) 


(AJ8) 


(AJ9) 


We  illustrate  g  and  in  Fig.  10. 


**)  {  x  (2/t*XQ/«^J 


-(«/»)!* 


s 


Figure  10.  The  pulse  M  of  Appendix  A.  Notice  that  g  Is  not  smooth  at  s  *  vl 2,  but  fit,*  is 


smooth  everywhere  in  (0,  *]. 


APPENDIX  B:  MATHEMATICAL  COMPLEMENTS  TO  SECTION  ID 
Proof  of  Theorem  3.2.  (i)  First  note  from  Eq.  (3.37)  that 

(dt^/dt.)#.)  -  1  -  [77G(ftt0)J[to  G(Gt0)]';  (B.l) 

since  Urn  [T/G(Qt„)]  «  ~  -  lim  [In  G(ftt0)]'  then 

l-*0* 

Um  (dtre/dto)(t0)  -  .  (B.2) 

Also,  sinrift  G/  and  G  are  continuous  on  (0,  x/fi)  and  G  never  vanishes  there  then  G'/G2,  hence 
is  also  continuous  there.  Therefore,  there  exists  8  >  0  such  that  (dt^/d^XtJ  <  0  for  all 
t0  €  (0, 5].  So  5  is  a  member  of  the  set  on  the  RHS  of  Eq.  (3.36)  and  t„  £  8  >  0,  as  required. 

(ii)  Since  IQ  >  0  then  [0,  TJ  *  0.  Let  to  e  [0,  TJ;  if  to  =  0  then  tD  is  admissible,  so 
suppose  tc  >  0.  If  C  <  t„  then  from  Fig.  2  we  see  that  0  #  P(t0)  c  I°(C ),  so  suppose  t  €  F(tc) 
and  note  that  dZ/dt  and  d^dt2  exist  at  (tj ;  t)  ami  (to;  t).  Now  from  Eq.  (3.9)  we  have 

(3/dt0)(32Z/9t2)(t0/ ;  t)  «  -(e/m.e^pG^C ;  t),tXdZ/dtJ(tJ ;  t)  (B.3) 

and  from  Eq.  (3.21)  and  the  fact  that  g(s)  *  0  only  for  at  most  s  =  0,  it,  it  follows  that,  for 

0  <  tj  <  t0,  LHS  of  Eq.  (B.3)  <  0;  but  p(Z(t^  ;t),t)  <  0  (since  charge  sheet  T0  is  at  Z(t„ ;  t)  at  time  t). 

Hence 


0Z«to)(to,;t)  <  0. 


(B.4) 


In  other  words,  Z(tD;  t)  is  a  strictly  decreasing  function  of  t0  at  each  t^  e  (0,  tD)  so  that  for  such  t/ , 
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Z0.';t)>Z(to;t)  . 


(B.5) 


But  also  Z(0;  t)  *  vct  >  Z(tG;  t)  so  Eq.  (B.5)  holds  as  well  for  t»'  *  0.  Hence  t0  is  admissible  by 
Definition  3.1.  Next,  if  T0  <  x/£l  then  there  exists  e  >  0  such  that  if  le  e  (t*  t„  +  e)  then  t0  is  not 
admissible;  indeed,  if  e  >  0  is  small  enough  (so  that,  at  least,  T0  +  e  <  x/Q)  it  is  clear  from  Fig.  2b  (or 
from  the  definition  of  TD  in  Eq.  (3.36))  that,  for  %  €  (T0,  T0+e),  sheet  te  leaves  after  and  returns  after 
sheet  T0  so  there  exists  t  e  F^J  n  I°(T0)  such  that  Z(l„;  t)  =  Z(T0;  t);  hence  ^  is  not  admissible. 
Finally,  if  To  =  nJQ,  then  [O.TJ  =  [0,  x/Q]  which  is  clearly  the  largest  possible  interval  of  t0’s.  ■ 

We  have  also  shown  (and  will  use  later  on): 

Corollary  B.1:  If  t0  e  [0,  fj  and  t  €  I°(t0)  then  the  map  t0  Z(tD;  t)  is  strictly  decreasing 

at  ta. 


We  next  prove  Theorem  3.3. 

Proof  of  Theorem  3.3.  If  (z,  t)  e  DJT,  G],  z,  t  >  0,  the  existence  of  a  solution,  say  tol,  in 
[0,  TJ  follows  immediately  from  the  definitions  of  D,[T,  G]  and  D[T,  G].  If  to2  is  another  solution 
in  [0,  tj  then  Zfl,,;  t)  =  z  =  Z^;  t).  But  t„j  and  t^  are  both  admissible  and  if  z  >  0  then  e  F(t0l) 
Fftoj);  hence  t*,  =  t^.  Next,  if  p(z,  t)  *0  then  there  exists  t^e  [0,  tj  such  that 
t)  =  zso  (ZfW  t),  t)  €  D,[T,  G],  i.e.,  (z,  t)  €  D,[T,  G];  hence  (z,  t)  e  D,[T,  G]  implies 
p(z,  t)  =  0.  Finally,  if  z  >  0  and  t  £  0  then  it  is  clear  that  p(z,  t)  =  0.  ■ 
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APPENDIX  C:  RETARDED  TIME 


In  tills  Appendix,  we  establish  the  existence  and  uniqueness  of  solutions  t ,(t0;  x,  t)  to  Eq. 
(4.58)  as  well  as  derive  smoothness  properties  of  these  solutions.  Although  our  results  may  appear  to 
be  physically  obvious,  we  eschew  proof  based  on  "obvious  physical  grounds".  We  fix 


*-($.0,0.  5>0.  C*0,  t£0,  and  t0  e  [O.TJ 


(C.1) 


throughout  unless  otherwise  noted.  (We  do  not  require  %  >  a  in  this  Appendix  C  only.) 
It  is  convenient  to  first  consider  the  case  t  =  0.  Eq.  (4.58)  then  becomes 


t  -t' 


d/c)K2  ♦  (v.t'-O2]1®,  if  t'  ;>  o 
l|x|/c,  if  t'  <  0 


(tB  =  0) 


and  tills  has  unique  solution 


(C.2) 


t'(0;x,t) 


K^Wv„H[t -(^((Ao)]2  -  H-(^)l][tJ-<^)J(|x|/v0)2l}w),  if  t  *  |x|/c 
t  -  |x|/c,  if  t  <|x|/c 


(C.3) 


where  the  first  line  may  be  rewritten  as 

"  (4>{[l-(i)J]^  +  (vjt-02}"*)  if  1 1>  |x|/c;  (C.3a) 

see  Fig.  11. 

We  next  consider  the  case  t0  *  0.  Rewriting  Eq.  (4.58)  in  detail  we  have 
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• 


/ 
: '  /*// c 


*'<P;  x,  t) 

ft)  i*  >  0  (f  > 


t) 


RgUre  1  **  7116  solution  to  Eq.  (C.2)  «riven 

^  *“«'%  uta.  In  {  “^'CaBy  ^  *»•  <c-3)-  R>r  UK  case  (b)  we 

•-Mb*  "■‘ItaeAfO 


Is  drawn  for  the 


(t.*0) 


(C.4) 


t-t' 


(m{?W -tj-ivjTM.at'W' -0K),}lfl.  if  t'  €  I(g 
|x j/c,  if  t;  g  i  (g 


where  we  have  used  Eq.  (4.71).  Unfortunately,  in  contrast  to  solving  Eq.  (C.2),  which  involves  as  an 
intermediary  only  a  quadratic  in  t',  solving  Eq.  (C.4)  involves  as  an  intermediary  a  quardc  in  t',  the 
solution  of  which  -  despite  the  availability  of  a  general  solution  -  is  extremely  tedious  to  obtain; 
furthermore,  we  would  still  have  to  sort  through  these  four  solutions  to  determine  which,  if  any,  are 
solutions  to  Eq.  (C.4).  (These  considerations  are  not  unwarranted.  For  example,  consider  solving 
(*):  aa-x  -  (a2-(x-a)2]1/2,  a  >  0,  a  e  K.  The  associated  quadratic  has  real  solutions  iff 
a  g  11-^2*.  1*^2].  Further,  if  a  g  [2,  then  both  of  them  solve  (*);  if  a  g  [0,  2)  then  only 
one  of  them  solves  (*);  and  if  a  g  [l-/f,0)  then  neither  solves  (*).)  Additionally,  this  t,,  *  0 
solution  -  if  Eq.  (C.3)  is  any  indicator  -  is  too  complex  to  be  integrated  analytically  in  Eqs.  (4.69) 
and  (4.70).  We  will  adopt  the  approach  of  not  solving  die  intermediary  quartic  but  ratter  of 
employing  geometric  reasoning  (as  in  the  solution  of  Eq.  (C.2))  to  establish  the  existence  and 
uniqueness  of  solutions  to  Eq.  (C.4);  flutter,  we  will  establish  the  smoothness  properties  that  we 
require  of  the  solution  without  ever  actually  obtaining  it  explicitly.  In  adopting  this  approach,  we  are 
abandoning  the  quest  for  expressions  explicit  in  the  inital  pulse  and  surface  parameters  for  the  finitely- 
ranote  Adds,  E,  and  B„.  We  are  willing  to  do  so  because  that  will  allow  us  to  further  progress 
towards  our  ultimate  goal,  namely,  obtaining  explicit  expressions  for  the  fields  in  the  radiation  limit, 
indeed,  it  win  turn  out  that  the  explicit  solution  for  V  when  t„  *  0  is  not  needed  to  obtain  the  explicit 
expressions  for  the  limit  fields. 

As  a  first  step,  we  graph  the  RHS  of  Eq.  (C.4);  the  curves  of  Fig.  12  include  all  generic 
possibilities.  The  curves  are  intentionally  drawn  to  be  "flat";  indeed,  denoting  the  function  on  the 
RHS  of  Eq.  (C.4)  by  JQt').  we  have,  for  tf  e  I°(t0), 


oT 


yo  uo 


Figure  12.  All  generic  possibilities  for  the  RHS  for  Eq.  (C.4).  Z*(to)  is  foe  maximum  z  value 
attained  by  sheet  tc  and  is  given  by  Eq.  (3.22)  as  ZM(t0)  =  (l/4)v077G(flt0). 
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Idj^dt'i  -  <ixi/c)|v0(t'-g  -  (v/rxxot^t ' - gJ  -  ci-|i-2(G(nt0yri(t'-g|/^(t/) 


<  ci)|i  -2[G<ot0)/7i(t/-  g|  <  i  <  1 


(C.5) 


since  |1  - 2[G(Qtoyr](t ' -  g|  <  1  for  t'- tc  e  (0,  77G(Q t„)) .  The  reader  should  note,  however, 
that  this  figure,  while  ^lightening,  will  not  be  the  basis  of  our  existence  and  uniqueness  arguments; 
rather,  they  will  be  independent  of  this  figure  and  based  solely  on  Eq.  (C.4). 

We  are  now  ready  to  prove  die  existence  remit 

Theorem  C.l :  Fix  x,  t  according  to  Eq.  (C.I)  and  tc  e  (0,  t#J.  Then  solutions  to  Eq.  (C.4) 

exist 

Proof.  First  define,  tor  if  £  t 

AGO-  t-t'  (C.6) 

and  suppose  first  that  0  £  t  <  |x|/c;  then  (Sf-  L,)(0)  >  0.  On  the  other  hand,  lim  L,(t')  »  ~  so 

*  l'-#— 

there  exists  t*  <  0  such  that  Lfi.")  >  |x|/c  «  SQX*);  hence  (if-  L)( t*)  <  0.  Since  if  -  L,  is 
continuous  on  (-«•,  °°)  then  there  exists  t  e  (t\  0)  such  that  (if-  L,Xt )  =  0,  so  that  t  <  t  solves 
Eq.  (C.2).  Similarly,  if  t  >  |x|/c  flien  (if  -  L,X0)  <  0;  and  L,(t)  ■  0,  t)  £  £/c  >  0  so 
(if  -  L,)(t)  >  0.  Hence  there  is  a  solution  in  (0,  t).  Finally  if  t  =  |x|/c  then  t/  =  0  <  t  is  a 
solution.  ■ 

While  existence  is  clear  geometrically,  simply  by  drawing  straight  lines  with  (any)  negative 
slopes  and  vertical  axis  intercept  (0,  t)  in  Fig.  12,  uniqueness,  on  the  other  hand,  is  a  bit  more  delicate. 
For  there  are  some  straight  lines  with  negative  slopes,  which  slopes  may  possibly  be  -1  since  we 
cannot  infer  much  from  file  figure  about  file  horizontal  (t')  and  vertical  scales,  that  intersect  the  curves 
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more  than  once;  and  for  such  lines  foe  solutions  to  Eq.  (C.4)  are  not  unique.  We  now  show  that  in 
fact  multiple  intersections  do  not  occur. 

Theorem  C.2:  Hxx,t  according  to  Eq.  (C.l).  Then  the  solution  to  Eq.  (C.4)  is  unique. 

Proof.  Suppose  that  L,  intersects  3[  at  two  distinct  points  whose  t7  coordinates  are  t/  <  t/ . 
Clearly  both  t/  and  t/  cannot  be  exterior  to  I(ta),  for  then  L,  would  have  slope  0.  Suppose 
t/ ,  t/  e  l(tc).  Since  ST  is  continuous  on  [t/ ,  il  ]  and  differentiable  on  (t/ ,  t/ )  then  by  the  Mean 
Value  Theorem  there  exists  t*  e  (t,\  t/)  c  I0^  such  that 

«LSj7dt'Xt*)  -  l^)  -  -  Ad/Mtj'-t/]  -  slope  of  if*  -l.  (C.7) 

But  we  have  already  seen  from  Eq.  (CJ)  that  |<t5f7dt y  |  <  1  on  I  “(tj  so  Eq.  (C.7)  yields  a 
contradiction.  Hence  one  of  t/ ,  t/  must  be  exterior  to  I(tD)  and  the  other  must  be  in  I(t„).  Suppose 
t/  <  t0  and  t/  e  I(tQ);  then,  since  dJ£7dt'  is  continuous  on  I°(t0)  (see  Eq.  (C.S)  and  note  that  iftt') 
is  never  0  since  £  >  0 ),  we  have 

«/ 

«f0a)  «  |x|/c  +  J(dJ^'dt/)dt/  (C.8) 


so 


(C.9) 


But 


|4,(t2')-£/t/)|  -  tZ-t/  >  t/-t0. 


(C.10) 
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Lfif)  -  |x|/c,  and  L(t/)  -  -SJ|( t,')  so  from  Eqs.  (C.9)  and  (C.10)  we  have 


M^t/)  -  {x|/c|  <  \S^)  -  |x|/c|  (C.11) 

which  again  is  a  contradiction.  Similarly,  supposing  t,'  €  Kg  and  tj  >  tjitj  also  leads  to  a 
contradiction.  Hence  L,  cannot  intersect  ^  at  two  distinct  points  and  the  uniqueness  proof  is 
complete.  ■ 


We  now  have  a  precise  enough  description  of  the  behavior  of  lines  L,(t  0  =  t  —  t /  relative  to 
the  curves  of  Hg.  12  to  use  this  figure  in  subsequent  proofs,  and  we  do  so  freely. 

The  next  result  in  this  section  will  be  crucial  in  evaluating  the  integrals  for  E,  and  B„  in  Eqs. 
(4.69)  and  (4.70),  giving  us  a  simple  description  of  die  set  of  all  tD  e  [0,  T„]  such  that,  for  given  x,  t, 
we  have  Xd  (*„;  t 7  (t„;  x,  t))  =  1,  thus  allowing  us  to  write  Eqs.  (4.73)  and  (4.74)  for  these  Adds.  For 
simplicity  of  notation  we  suppress  die  x  and  t  dependence  of  i'  and  write  simply  1 7 (tD).  Also,  we 
introduce  the  following  notation  for  die  various  boundaries  of  D[T,  G]: 

a»D[r,G]  -  (J  {(go), 

«.«tO.TJ 

a^ir.Gj-  (J  (Mjg)}, 

dLD[r,G]  -  {0}  x  [0,oo), 

and 

djDIT.G]  -  {T0}  x  I(g;  (C.15) 


(C.12) 

(C.13) 

(C.14) 
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additionally,  for  future  use  we  denote 


a^ir,  G]  «  a,D[r,  G]  u  a,D[r,  gi  .  (c.  i6) 

Theorem  C.3:  Hx  x,  t  according  to  Eq.  (C.l). 

(i)  If  -«»  <  t  -  |x|/c  S  0  that  for  all  ta  e  [0,  Tj,  t /(t0)  =  t  -  |x|/c;  so 

(t».  t ' (t„))  €  D ~[T,  G]  unless  t  -  |x|/c  =  0  =  t0  in  which  case  (t0,  t 7 (tc))  e  3„D[ T,  G]. 

(ii)  If  0  <  t  -  |x|/c  ^  Tc  then  (t  -  |xj/c,  t  -  |x|/c)  e  ^Dfr,  G], 

(tc.  t'(t„))  €  D°[r,  G]  whenever  tn  e  (0,  t  -  |x|/c),  and  (t0,  t'(tc))  e  D'|T,  G]  whenever 
t0  €  (t  -  |x|/c,  T0].  (If  tD  =  t  -  |x|/c  then  t'(t0)  =  t  -  |x|/c  so  (t0,  t'(tc))  e  3,(7,  G].) 

(iii)  If  Tc  <  t  -  |x|/c  <  tn.CTj  then  (Tc,  t  -  |x|/c)  €  3*  D [T,  G]  and 

(to.  t'(t0))  e  D°[r.  G]  whenever  t0  €  (0, 10);  also  (I0,  t '(!„))  e  {Tj  x  r(T0)  C  D[7,  G]  but  it  is 
never  true  that  t'(to)  =  t  -  |x|/c. 

(iv)  If  t^IJ  St-  |x|/c  <  «  then  (t,(t-|x|/c),  tj(t-|x|/c)]  x  {t— |x|/c>  c  Sj-Dfr,  G]  where 

t,(t-|x|/c)  ■  inf  C[t-|x|/c],  (C.17) 

ts(t-|xj/c)  *  sup  C[t-|x|/c},  (C.18) 


and 


t«B [t~|x|/c]  ■  {t0  e  (0,g|tflB(to)  -  t-|x|/c};  (C.19) 

also  Clt-lxl/c]  0  and  tKt  -  |x|/c)  >  0.  In  addition,  (t0,  t '(t0))  €  D°[T,  G]  whenever 
tc  e  (0,  -  |x|/c))  and  (tc,  l'(tc))  e  D17,  G]  whenever  tc  e  (ts(t  -  |x|/c),  T0]. 
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(If  t0  e  [t,  (t-|x|/c),  tj  (t—  jx |/c)J  then  t '(t0)  =  t  -  Jx|/c  so  (t0,  t ' (t„))  e  df  D[T,  G].) 

(v)  In  cases  (U)  -  (iv).  (0.  t '  (0»  e  dL  D(r.  G]. 

Proof,  (i)  It  is  clear  from  Figs.  11  and  12  (die  latter  augmented  with  the  line 
Ltf) « t  - 1*)  that  t/(t0)  =  t  -|x|/c  whenever  t  -  |x|/c  £  0  and  t0  £  0;  further,  unless 
t  -  |x|/c  *  0  «  t0,  then  t '  (te)  <  0  so  from  Fig.  2  we  have  (tc,  t '  (t„))  e  D'[r,  G],  while  if 

t  -  |x|/c  =  0  =  t0  then  (tc,  t'(t0))  =  (0,0)  e  a.Dfr,  G]. 

(ii)  Let  0  <  t  -  |x|/c  £  T0;  then  from  Fig.  2  we  have  (t  -  |x|/c,  t  -  |x|/c)  e  BgDir,  G], 

Let  t0  e  (0,  t  -  |x|/c),  i.e.,  t  -  |x|/c  >  tc;  then  supposing  t'(t„)  £  tc  gives 


Lfi'  Cg)  =  t-t'(t„)2>t- 10>  |x|/c  (C.20) 

while,  from  Eq.  (C.4)  or  Fig.  12,  *  |x|/c  so  that  f 7 (t0)  cannot  be  die  retarded  time  for  t0> 

which  is  a  contradiction;  and  supposing  t /(t0)  £  *  t0+77G(Qt0)  gives 

Lfi'w  -  t -l'<g  £t-[to  ♦  77G(Qg)  <  t - [f0  ♦  77G(flf0))<  t -T„  <;  |x|/c 

(C.21) 

while  again  *  |x|/c,  a  contradiction.  So  l7(g  e  (tD,  t^g)  *  I°(g  and  so 

(t„,  t7(t0))  e  D°[r,  G].  (Alternatively,  one  may  argue  using  Fig.  12  that,  when  0  <  t„  <  t  -  |x|/c  £ 

T0,  the  line  Ltf)  =  t  - 1 '  intersects  any  of  die  curves  at  some  tf  e  (tc,  TJ  c  (t„,  i»(t0»,  since  T0  < 
tBB(t0)  for  all  tG  e  [0,  TJ,  so  that  t'(tD)  e  (tc,  (tc)).)  Now  let  t„  €  (t  -  |x|/c,  T0],  i.e.,  t  - 
|x|/c  <  t„;  then  it  is  clear  from  Fig.  12  that  t 7  (t0)  =  t  -  |x|/c  so  (t0,  t '(g)  €  D"[F,  G]. 

(iii)  Let  T0  <  t  -  Jx|/c  <  ^(T,,);  then  from  Fig.  2  we  have 
(T0,  t  -  }x|/c)  e  dnD[r,  G].  Let  tc  e  (0,  T0);  then  supposing  t 7  (t0)  <,  te  gives,  similarly  to  Eq. 

(C.20), 


98 


*  t  -t„  >  t -T0  >  |x|/c  -  ^tf'(g). 


(C.22) 


•  contradiction;  and  supposing  f '(t0)  S  gives,  similarly  to  Eq.  (C.21), 

Lfi'itJ)  <  <  t  - (t  -  |x|/c)  -  |x|/c  -  '(g),  (C.23) 

also  a  contradiction.  So  (tc,  t'(t0))  e  D°[7*,  G].  Also,  that  (Tc,  t / (T0))  e  {T0}  x  F(T0)  may  be  seen 
from  Fig.  12  by  taking  t0  =  T0  there;  note  also  that  t  '(T0)  <  t  -  |x|/c  or  I ' (t0)  >  t  -  |x|/c  so  that 
it  is  never  true  that  t  ;(t0)  =  t  -  )x |/c. 

(iv)  Let  g  (I0)  <  t  -  |x|/c  <  «;  then  from  Fig.  2,  allowing  for  the  possibility  that  t^g  has 
subintervals  of  constancy  in  (0,  TJ,  we  have  [t,,  y  x  (t  -  (x  |/c )  q  djD[T,  G]  where  we  have 
suppressed  the  argument  t  -  |x|/c  of  t,  and  y  Since  t^tj  increases  from  t^CT,,)  to  <»  as  t0 
decreases  from  T„  to  0,  then  t^[t  -  |xj/c]  *  0  for  t.JI,,)  St-  |x|/c  <  ~;  so  ^  >  0.  Let 
t„  e  (0,  t,);  then  supposing  t '  ( to )  S  t0  gives,  similarly  to  Eq.  (C.20), 

Lfi'W  £  t  -tQ  £  t  -T„  >  t  -tjfj  ;>  |x|/c  =  ^a'(g),  (C.24) 


a  contradiction;  and  supposing  t'(t0)  2:  tj(t0)  and  denoting  by  t*  any  member  of  [t,,  y,  so  that 
trtn(to)  =  t  -  |x|/c,  gives 

A(t '(O)  <  t  -  tjy  <  t  -tjt,)  S  t  -tm(0  =  t  -(t-|x|/c)  =  |x(/c  =  (C.25) 

a  contradiction.  So  (tc,  t'(t0))  e  D°[r,  G].  Now  let  t0  e  (ts,  TJ;  then  with  t*  as  above, 

^(Uto)) -t -ug  >  t -ug  ;>  t  -uo  - 1  -(t-|x|/c)  =  |x|/c  (c.26) 
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while  LJt)  =  0.  But,  from  Eq.  (C.4)  we  have 


^tjg)  -  |x|/c  -  fit),  (C.27) 

the  second  equality  holding  because 

t  >  t  -  |x|/c  «  tjO  k  tm(g  >  tjfj.  (C.28) 

So  Litf)  and  if(t ')  must  have  their  unique  intersection,  which  occurs  for  t7  *  t '  (t„),  at  a  value  of 
t'  €  (t,*  (ta),  t);  i.e.,  t'(t0)  >  Ut„)  and  so  (tc,  t '  (t0))  e  D-[r,  G). 

(v)  Follows  from  Fig.  1 1  by  noting  that  I '  (t„)  >  0  whenever  t  -  |x  |/c  >  0.  ■ 

The  results  are  represented  diagrammatically  in  Fig.  S,  which  we  referred  to  earlier  in  Section 
IV. 

The  above  mathematical  result  has,  of  course,  physical  interpretation.  For  example,  (i)  says 
that  if  die  required  retarded  time,  namely  t  -  |x|/c,  for  a  charge  sheet  located  at  z  =  0  is  negative 
then  no  sheet  can  satisfy  that  requirement;  this  is  of  course  clear  since  no  sheet  is  in  flight  at  any 
negative  time.  The  other  parts  have  similar  interpretations;  we  leave  their  explicit  elucidation  to  the 
reader. 

We  remark  here  that  the  accommodation  in  the  last  theorem  of  the  possibility  of  intervals  of 
constancy  of  t^  is  not  unwarranted  generalization  --  we  presented  in  Appendix  A  a  class  of  pulses  for 
which  tm  has  such  an  interval  of  constancy. 

Theorem  C.3  tells  us  that  the  map  tD  **  Xo  (t0;  t /  ( t0 ;  x,  t)),  for  fixed  x,  t,  is  piecewise 
continuous  on  [0,  T0],  being  given  by 
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X0(to;t'(to;x,t)) 


1,  if  t„  e  [0,To(x,t)] 
0,  if  tc  €  (T0(x,  t),  TJ  , 


where  T0(x,  t)  is  given  by 


T0(x,  t)  ■ 


0, 

t-|x|/c, 

V 

t,(t-|x|/c). 


if  -«»  <  t  -  |x  |/c  £  0 
if  0  <  t-|x|/c^T0 
if  T0  <  t  -  |x|/c  <  tJTj 
if  tJO  *  t-|*|/c  <  «  . 


(C.29) 


(C.30) 


Hence  the  %x>  factors  in  Eqs.  (4.69)  and  (4.70)  can  be  accounted  for  simply  by  choosing  0  for  the 
integration  lower  limit  and  T0(x,  t)  for  the  upper  limit  and  setting  Xd  *  1-  In  Fig-  6  we  plot  T#(x,  t) 
as  a  function  of  t  for  fixed  x,  with  one  jump  corresponding  to  one  Interval  of  constancy  for  t^tj,  as 
in  Fig.  5;  recall  that  we  also  referred  to  Fig.  6  earlier,  in  Section  IV. 

Strictly  speaking,  the  upper  limit  for  foe  last  case  of  Eq.  (C.30)  should  be  t$(t  -  |x|/c)  (see 
Theorem  C.3(iv».  We  have  chosen  rather  than  tj  because  the  choice  of  t,  allows  a  more  unified 
(Le.,  caseless)  treatment  in  foe  sequel  than  does  foe  choice  of  V  This  point  is  almost  always  moot 
since  in  almost  all  non-contrived  cases  we  have  tj  =  tj  (the  pulses  of  Appendix  A  are  examples  of 
"contrived"  pulses).  Nevertheless,  we  will  be  careful  to  point  out  and  include  those  modifications  to 
our  formulation  that  are  necessary  in  case  t*  >  t,.  As  a  point  of  interest,  note  that  foe  choice  of  t$ 
rather  than  t,  for  the  definition  of  Tc  would  make  the  map  t  >-»  T9  (x,  t)  left-continuous  rattier  than 
right-continuous  at  the  jump  in  Fig.  6. 

The  last  two  results  of  this  Appendix  establish  the  continuity  and  differentiability  properties  of 
the  map  t„  ~  t^t,;  x,  t),  with  x,  t  fixed  and  tD  variable  according  to  Eq.  (C.l),  that  we  will  need  in 
Appendix  D.  Consistent  with  our  aforementioned  convention  of  suppressing  x  and  t  in  the  argument 
list  of  ,  we  also  write  (dt  Vdt^O  to  mean  (9t//9t0Kto;x,t).  We  denote 
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D t[T,  G)  ■  [0,  TJ  x  ~), 


(C.31) 


with  interior  given  by 


D:[ r,  G]  =  (0.  Tc)  x  (-«o,  °°);  (C.32) 

and  we  also  denote 

(D^jr,  G]  ■  D^r,  G]  \  a^Dtr,  gi  -  D°ir,  G]  u  d -°[t,  gj  (c.33) 


(see  Eq.  (C.16)). 


Theorem  C.4:  Fix  x,  t  according  to  Eq.  (C.l). 

(i)  If  -«  <  t  -  |x|/c  S  0  then  t /(t0)  is  continuously  differentiable  (C1)  on  [0,  TJ. 

(ii)  If  0  <  t  -  |x |/c  S  fc  then  t,(t0)  is  Cv  on  (0,  t  -  |x|/c)  u  (t  -  |x|/c,  T„)  . 

(iii)  If  T„  <  t  -  |x|/c  <  U(To  )  then  t(g  is  Cl  on  (0,  T„)  . 

(iv)  If  t^CT^  St-  |x|/c  <  <*>  then  is  C  on 
(0,t,(t-|x|/c))  u  (t,(t—  |x |/c),  tj(t”|x |/c))  u  (ts(t— (x |/c), T„) . 

Further,  for  those  tQ  at  which  t /{t0)  is  C1  we  have  in  each  case 


(dt'/dt^g 


- ( l/cJACg  (d^dtoKt,,) ,  if  t0  €  FI  [a  first  interval  of  (i)-(iv)] 

0,  if  t0  e  FI  [not  a  first  interval  of  (J)-<iv)  ] 


and 


(C.34) 
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(dt'/dg(g 


Ki/c^goz/ago.;  *  '(UM  i  •Ki/c)A(t(>xaz/dtxt0;  t  '(g)], 

o. 


if  t„  6  FI 


ifto€  FI 


(C.35) 


where 


8(g  -  z^f'cg). 


(C.36) 


or,  more  precisely. 


8(t0;  x,  t)  *  Z(t0;  t  '(t0;  x,  t)) , 


(C.  37) 


and 


A<g  « [8(g  -  C)/{?  ♦  I8(g-C]3}w  (c.38) 

(and  the  arguments  of  Z  are  denoted  (t0,  x)). 

Proof.  The  proof  of  (i)  here  follows  immediately  from  (i)  of  Theorem  C.3,  for  we  have 
dt7/dt0  *  0  on  [0,  TJ.  To  prove  (ii)-(iv)  we  appeal  to  the  Implicit  Function  Theorem  which  says  that 
if  real-valued  function  TG(t0,  xf)  has  continuous  first  partial  derivatives  dtH/dr,  and  dM/dxf  in  some  open 
subset  of  R1  which  contains  interior  point  (t0,  tj )  satisfying  0-C(to,  tj )  =  0  and  (dTCGt'Xt,,,  t^ )  *  0, 
then  there  is  some  nonvoid  neighborhood  (t0  -  5,  to  +  5)  of  to  and  a  unique  function  h  defined  on 
(t„  -  5,  t0  +  8)  such  that  t0'  =  h(t0)  and  Oifo.  h(x„))  =  0  whenever  t„  e  (t0  -  5,  t0  +  6);  further,  h  is 
C*  on  (t„  -  5,  t„  +  5).  To  apply  this  to  our  present  situation  we  first  define,  for 
(t„,  t')  e  D:tr,  G], 
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tK(t0, t ')  -  t'-  t  +  J£t'). 


(C.39) 


We  next  note,  using  a  rewritten  version  of  the  RHS  of  Eq.  (C.4),  that 


(l/c){§2  +lZ(t0;  1 0  -CP } w. 
|x|/c. 


if  t '  e  I(t„) 

if  t'  6  I-(g 


(C.40) 


and,  from  the  differentiability  properties  of  Z(t0;  t)  on  D([7,  G]  (see  Eq.  (3.34)),  that  d90dto  and 
tfHfdl'  exist  and  are  continuous  on  (D“)k[7,  G]  with 


0W3t,Xto.t/) 


1  ♦  (d^dt'Xt').  if  t'  e  Itg 
i,  if  t'  e  r(g 


(C.41) 


there;  and  we  then  use  Eq.  (CJ)  to  conclude  that  Wfd\‘  *  0  there.  Now  if  to  e  (0,  To)  then  our 
previous  existence  and  uniqueness  results  imply  flat  (t„,  f  ,(tD))  e  D“[7,  G]  is  such  flat 


M(t0,t,(t0))  =  Oforallt0€  (0,To) 


(C.42) 


(not  just  tD -locally)  so  must  be  the  unique  h  guaranteed  by  the  Implicit  Function  Theorem 
whenever 

M'<g)  €  Ws  [T’G]  ;  (C.43) 

hence  I7  must  be  C1  at  such  t0.  But  items  (iiHiv)  of  Theorem  C.3  indicate  that  Eq.  (C.43)  is 
satisfied  precisely  for  those  t0  specified  in  cases  (iiXiii)  here,  and,  in  case  (iv)  here,  for 
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t0  e  (0,  tj  u  (4,  T„).  Since,  in  case  (iv),  t /  is  in  fact  constant  on  [t,,  4]  (being  t  -  |x  j/c  there)  then 
we  may  include  the  interval  (t,,  4)  so  that  (iv)  is  also  proved.  Lastly,  the  first  expression  far  dtVd^ 
results  from  differentiating 


Hgf'O.))  -  t'(g  -  t  +  Jflf '<g) 


(C.44) 


(see  Eq.  (C.39))  with  respect  to  t0,  using  Eq.  (C.42);  and  the  second  expression  for  6t' /dto  results 
from  noting  that 

(d3/dg(g  -  @Z/dg(to;t'(g)  ♦  0Z/&t)(to;t/(to)X(ff,/dto)(g.  (C.45) 

In  both  cases,  the  continuous  differentiability  of  Z  on  (DJDJT,  G]  and  of  t/  on  the  intervals 
in  (IMiv)  guarantee  die  existence  of  (continuous)  d8/dt0  on  those  mw  intervals  ■ 

Hie  above  theorem  indicates  that  there  may  be  exceptional  points  in  [0,  Tj  where  fails  to 
be  C*.  The  following  is  true  however. 

Theorem  C.5 :  Fix  x,  t  according  to  Eq.  (C.l).  Then  t '  is  continuous  on  [0,ToJ. 

Proof.  If  -°°  <  t  -  |x|/c  0  then  file  result  follows  immediately  from  Theorem  C.4(i).  So 
let  0  <  t  -  jx|/c  <  ~  and  consider  first  continuity  at  t0  =  0. 

Define 


C  »  min  {t-|x|/c,  t,(t-|x|/c),T0) ; 


(C.46) 
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by  items  (UMiv)  of  Theorem  C.3,  it  follows  that  t '  (t„)  e  P(te)  for  all  to  e  (0,  C)  and  so  for  such 
to,  jgj '(g)  is  given  by  the  first  line  of  Eq.  (C.40).  We  claim  -  and  show  below  -  that  lim  f '(t0) 
exists;  then  lim  Jf^t  /(t0))  exists  and  is  given  by 

lim  JRl'O  -  d/c){^  ♦  K  lim  t'(g  -  d2}"2  ,  (C.47) 


since  lim  G(Qg  ■  0  and  the  square  and  square  root  are  continuous  on  non-negative  reals.  Hence 

t-»o*  * 

from  Eq.  (C.44)  we  have,  using  (from  Eq.  (C.42))  lim  tK(tc,  t/(t0))  =  0,  that 


t  -  lim  f '<g  -  (l/c){$2  «•  [vc  lim  t'(g  -  CHW.  (C.48) 

v— o* 

where  t  ,(tc)  >  t0  2  0.  But  this  last  equation  is  foe  same  as  that  of  (the  first  line  of  )  Eq.  (C.2)  and 
we  know  the  latter  has  unique  solution  I '  (0)  given  by  Eq.  (C.3).  So  we  must  have 


lim  t'(g  -  f '(0),  (C.49) 

v-o* 

i.e.,  I/(t0)  is  continuous  at  t0  =  0. 

To  show  that  lim  t7(t0)  indeed  exists  we  proceed  as  follows,  still  requiring  t0  e  (0,  C).  In 
Eq.  (CAS)  for  d3AJt0  we  substitute  the  first  line  of  Eq.  (C.34)  for  dt  7/dt0  to  get 


r 
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(d3/dt.Xto)  “  0Z/dto)(to;t/(to))  /  [l  ♦(l/c)A(t0XdZ/<hXt0; ! /(t0))] 


(C.50) 


where  we  can  be  sure  that  the  denominator  on  the  RHS  is  strictly  positive  since,  by  Eqs.  (C.38)  and 
(3.35), 

|d<g|  <  1  and  |(dZ0t)(to;t/<to))|  £  v#.  (C.51) 

But  (3Z/9to)(to; t)  <  0  for  (to,  t)  €  D°[r,G],  as  per  Corollary  B.l,  so  that  <  0; 

hence 

(dS/dt^g  <0,  tc  €  (0,0.  (C.52) 

i£.,  8  is  strictly  decreasing  on  (0,  C).  And  by  Eq.  (C.34), 

(dt'/agoj  -  (i/c)|(d8/dt0xg|A(g,  tD  e  (o,  O-  (C.53) 

Now  if  s  m  sup  JKg  £  C  then  A(g  £  0  on  (0,0  so  that 

(dt'/dg(g  £  0,  5  £  C,  tc  e  (0,0  (C.54) 

while  if  $  >  £  then,  because  3  is  strictly  decreasing  on  (0,0.  there  exists  C*  e  (0,0  such  that 
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A(t„)  >  0  on  (0,  C)  *o  that 


(<€/AHcKt0)  >0,  s  >  C,  tc  e  (0,C).  (C.55) 


Further,  since  0  <  t„  <  t  '(t0)  <  t  on  (0,  C),  then  in  either  esse  t  is  bounded  and  monotonic  on 
(0,0  so  that  Um  exists  and  our  claim  is  proved. 

Lastly,  we  show  continuity  at  the  other  exceptional  points.  To  that  end,  let  be  defined  by 


t  -  |x|/c, 

t,(t-|x|/c). 


if  0  <  t  -  jx|/c  £  T0 
ifT„<t-|x|/c<tIIB(f0) 
if  *  t-jx(/c  < 


(C.S6) 


then,  for  t0  e  (0. 0>  Eqs.  (C SI)  and  (CJ3)  still  hold.  Now  if  i  ■  inf  8CU  £  C  then 

«.«(aO 

MtJ  2;  0  on  (0, 0  so  that 

(df'/tkxy  *  o,  i  ;  tc €  (0, o  (C.57) 

while  if  i  <  5  then  ffiere  exists  C*  €  (0, O  such  that  A(tc)  <  0  on  (O  O  so  that 

(tf'/dt^y  <0,  i  <  C,  tc  o  CO-  (C.58) 
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Hence  i'(t0)  is  bounded  and  monotonlc  on  (C*.  O  so  lim  t  (t„)  exists.  Then  from  Eq.  (C.44)  we 
have 


t  -  lim  f'(g  -  (l/cjfc+fvj  lim  l ' 

V*<0*  (  «.-0* 


(g  -  ci  -  iv/rwaoi  liming  -ci 

«.-K Cr 


wj1 


(CJ9) 


so  by  of  solution  to  this  equation  (which  is  the  first  line  of  Eq.  (C.4))  we  have 


Umt'cg  -l'(0.  (C.60) 

'.-♦(O' 


This  is  the  continuity  result  in  case  C  *  T„;  otherwise,  we  have  from  Theorem  C.3  that 


lim  !  /(to)  ■  t  -  |x|/c  «  l '(C)  (C.61) 

«.-<0* 
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lim  l'(g  -  f'(0  -  t  -  jx|/c  (C.62) 

v< 


which  is  the  continuity  result  for  the  other  two  values  of  C  (as  wen  as  for  C  *  t*(t -  |x|/c)). 
HnaUy,  continuity  at  t0  *  Tc  in  cases  t  -  |x|/c  £  Io  and  t  -  |x|/c  2  tnn(Io)  follows  from  the  fact 
that  f '(to)  «t-  |x|/c  in  some  (one-sided)  non-punctured  neighborhood  of  T0.  ■ 
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Corollary  C.6:  Fix  x.  t  according  to  Eq.  (C.l).  Then  the  map  to~  8(to;  x,  t)  is  strictly 
decreasing  on  (0,  T0(x,  t)]. 

Proof.  From  Eq.  (CJ2)ff,  d8/dtc  <  0  on  (0,  T„(x,  t));  and  from  Eq.  (C.37)  and  Theorem  C.5, 
the  map  is  continuous  on  [0,  T0(x,  t)].  ■ 


HO 


APPENDIX  D:  PROOF  OF  THEOREM  5.1 


We  present  our  work  as  a  sequence  of  five  results  culminating  in  existence  Theorem  5.1.  We 
will  use  in  the  sequel  A.  rewritten  in  the  form 


A.(tc;t)  -  (1/sin*)  |x|-‘  Z(t0;t)  -  cot*. 


(0.1) 


Lemma  D.l :  Let  a,  b,  c„  c*  d„  dj,  d  e  ft,  a  <  b,  c,  <  dj  (i  »  1,  2),  c'  >  0;  and 
M,  ■  max{  |c,|,  |dj|}  >  0.  Let  functions  h :  [a,  b]  x  [Cj,  d,]  ->  ft  and  [a,  b]  x  [d,  ~)  -* 


(Cj»  dj  be  such  dial 

0)  h  is  continuous,  with  lAl„»M<«*; 

(ii)  Vy  €  (c,  m)  (ta;  y)  -  *,(t0;  y)  is  continuous,  i  =  1, 2;  and 
(Hi)  V,  €  M  lim  *,(t0;y)  -  ^.(t^  e  [c,,d,J,  i  -  1,  2,  uniformly  on  [a,  b], 

•  y-»- 

Further,  for  y  e  [a,  b]  define 

b 

3%)  -  J  <»„  *j(t0;y)W0,*,(t0;y)). 

■ 

Then 


o 

lim  !K(y)  «  fdt,  t„, 

y~*~  i 


(P2) 

(0.3) 


Proof.  Let  e  >  0.  Since  h  is  uniformly  continuous  on  [a,  b]  x  [c„  d,J,  there  exists 
8>  0  such  that  Vl€fe„  |*  -*'1  <  8 ->  |A(t„,*)  -  fc(t0,*/)|  <  e/2Mj(b-a).  Since  the 
convergence  of  *,(t»;  y)  to  ^.(U  is  uniform,  then  3Y>0  V,  €  feb,  y  >  Y  ->  |*,(t0;  y)  -  *,  _(t„)  |  <  8 
and  |kj(t0; y)  -  *l_(t0)|  <  e/2(M  +  l)(b  -  a);  further,  is  continuous,  hence  integrable,  on  [a,  b]. 


Ill 


Now  if  y  >  Y  then  V,  . 

|*J..(to)h(to,*1..(t0))  -  *I(t0;y)/«<t0,*1(t0;y))| 

*  l*(t0.Mto))l'l*2..(t0)  -  ¥{o'  y)  I  -  IVto:y)H*(to**..-(to)>  - 

£  (M  ♦  l)*e/2(M  ♦  l)(b  -  a)  ♦  Mj-e/2Mj(b  -  a)  =  e/(b  -  a) 

and  so 


b 

-  k2(te;yMe,ti(t0;y))l  <  e 

■ 

whenever  y  >  Y.  ■ 

Proposition  D.2 :  Let  *  -  (sin$,  0,  cos$)  with  0  £  $  <,  n 12 ,  T  >  0,  and  tD  g 

Then 


lim  |* I*1  Z(t0;f'(t0;|xj*,  |xf/c  «•  t))  «  0 

Vuniformly  on  [0,T“(t)]. 

Proof.  From  the  first  line  of  Equation  (C.3),  which  applies  here  since  jx|/c  + 
compute 


Um  l'  (0;  |x|&,  |x|/c  ♦  x)  *  x/[l  -  (v  /c)  cos$]. 
M— 


(D.4) 


(D.5) 


[O.f „•(!)]. 


(D.6) 


x  >  |x|/c,  we 


(D.7) 
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Since 

ZCO-.t'cO;  |x|£.  |x|/c  ♦  x))  «  v/(0;  |x|£.  |x|/c  ♦  x) 

then 

Um  |*rZ(0;  t^O;  |x|£,  |x|/c  +x))  -  0. 

1*1-*- 

But  by  Corollary  C.6, 

0  S  Z (t„; t (tG;  |x|£,  |x j/c  ♦  x))  <;  Z(0;  t'(0;  |x|£,  |x|/c  +  x)) 

for  all  to  e  [0,To"(x)]t  since  T0(|x|£,  |x|/c  ♦  x)  «  T”(x).  Hence,  we  arc  done.  ■ 

Proposition  D.3:  Let  £  «  (sin$,  0,  cos$)  with  0  <  $  S  n/2,  x  >  0,  and  tc 

Then 

(i)  we  haw 

Um  t'cg  |x|£,  |x|/c  +  X)  =  tl(t0;<^,x) 

1*1-*- 

Vpointwise  on  [0,  T”(x)] ,  where 


(D.8) 


(D.9) 


(D.10) 


[0,T"(x)]. 


(D.ll) 
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x/C,  (♦). 


if  t  -  0 


t»(t0;4>.x)  - 1 


x. 


x  -  [l/C2(<fr)G(QQ] ((:,(*) 

♦  C2(^)G(Qt0)(x  -  g 
-  icfo)  +  2C2(<»G(G0(t  -  gi1/2}. 


if  tD  *  0  and  $  =  i02 


if  t0  *  0  and  0  <  $  <  n/2 


(D.12) 


for 


C,($)  -  1  -  (v^c)  cos$  and  C2($)  -  27’*1(v</c)cos$ ;  (D.13) 

(ii)  to  -» t'jg^.x)  is  continuous  on  [0,T“(x)]; 

(iii)  the  convergence  in  (i)  is  actually  Vuniform  on  [0,  T  “(x)] ;  and 

(iv)  lim  tl(t„;^,x)  =  x  for  all  t0  e  [0,f  "(x)]. 

♦-K«ar 

Proof.  (i)  If  ^  =  0,  dirai  the  result  follows  immediately  from  Equations  (D.7)  and  (D.13).  so 
suppose  t,  >  0.  Since  x  >  0,  then  Equation  (S.16)  holds  for  all  |x|  >  0  so  that,  from  Equation  (C.4) 
with  x  *  |x|&  and  t  *  |x|/c  +  x,  we  have  for  all  |x|  >0 

|x|/c  «•  x  -  f'(t0;|x|*,  |x|/c  +  x) 

(D.14) 

■  (|x|/c){sin^  +  [|x|-‘Z(to;t/(to;  |x|«,  |x(/c  ♦  x))  -  cos2*]2}”2. 

Squaring  both  sides  of  this  equation  and  using  t'-t,,  *  (t ' -%)  *  (x  -g  yields 
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CJWG(Qt0)I!/(t0;|x|i,  |x|/c  +  x)  -x]2 

♦  2[C,($)  ♦  Cl(<»G<ag< x  -  to)J[t/(t0;  jx Jft,  |x|/c-x)  -  t] 

-  (t  -g[2  -  2C,«»  -  C2(^)G(Qt0Xx  -g]  (D.15) 

-  D„(t0;|x|)  +  DjCglxDlfcglxIfi,  |x|/c  +  x)-x]J 

♦  D^glxD^glxli,  |x|/c  ♦  x)  -  Xl  ♦  D0(to;  |x|) 


where  C,(4)  and  C2($)  are  given  by  Eq.  (D.13)  and 

DM<gl*l)  *  -<|x|/c)*'(v</c)17,*1G(flg[f/(t0;|x|«,  |x|/c  +  x)  -  x]3 
x  {7,*,G(Qt0)|t/(to;  |x|&,  |x|/c  +  x)  -  xj  +  2[2T-‘G(Qt0Kx-g-l]}, 


(D.16) 


D2(to;|x|)  -  (|x|/c)“(l  - (V</c)J  {1  -6T-1G(nt0Xx-g  ♦  flT-^QgCx-gj*}),  (D.17) 
D,(t0;  1*1)  -  -2(|x|/c)-|(vf/c)2(x -g{i  -sr-^ngfx -g  ♦  2[r-,G<at0xx -gi2},  (D.i8) 


and 


D0(to;  |x|)  •  -i\x\/cyl(vjc)\z  -gJn  -r-KKogfr-gj2 ;  (D.i9) 

note  that  we  have  supressed  the  x-dependence  of  Dq,  Dj,  D2,  and  D„  and  the  $  -dependence  of  DM  in 
their  argument  lists.  In  the  above  we  have  introduced  extraneous  roots  for  f  in  the  squaring  process 
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but,  by  Theorems  C.l  and  C.2,  Eq.  (D.14)  has  unique  solution  for  all  |x|  >  0  and  t  >  0.  Now 
since  tc  e  10,  f0"(t)]  then  t  €  [t0,  t0  +  T/G(£J  g] ;  combining  this  with  Eq.  (5.16)  we  then  have  for 
all  |x|  >  0 

|f/(t0;  |x|i,  )x |/c  +  x)  -t  |  <  T/G(ftte).  (D.20) 

Thus 


|D«(t0;  |x  |)  |  7(|x|/c)-‘  (V</c)2  T2/G2(Qg,  (Dll) 

|D,(t0; |x|)|  •|t/(t0;|x|«,  |x|/c  ♦  x)-x|2  £  (|x|/c)-*[l  ♦13(v,/c)2]r2/G2(Qt0),  (Dll) 

|D,(t0;|x|)|-|f(t0;|x|«,  |x|/c  ♦x)-x|  S  ndxl/O-^vyc^rVG^ftg,  (D.23) 

and 

|D0(to;  |x |) |  £  4(|x|/c)-1(vt/c)2ir2/G2(nt0)  (D1A) 


so  that 


|RHS  of  Eq.  (D.15)|  £  [1  ♦36(v,/c)2](|x|/c)-,r2/G2(Qt0)  <  2(|x|/c)'1r2/G2(flt0).  (D.25) 


If  $  -  */2  then  from  Eqs.  (D.15)  and  (D.25)  we  have 
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It'CV,  |x|t,  |x|/c  ♦  x)  -x  I  «S  (|x|/c)-‘  rJ/GJ(Qt0) 


(D.26) 


so 

limt'(t0;  |x|£,  |x|/c+x)  *  t  ($  -  x/2)  (D.27) 

Vpointwise  on  (0,  T"(x)] .  Since  this  result  must  hold  for  all  roots  of  Eq.  (D.15),  it  also  holds  in 
particular  for  the  unique  root  of  Eq.  (D.14),  so  the  result  stated  for  $  *  n/2  in  Eq.  (D.12)  follows. 

Suppose  then  that  0  <  4  <  x/2  and  consider  the  quadratic  equation  LHS  of  Eq.  (D.15)  =  0,  i.e„ 

CjWGfogw2  ♦  2[c,(4)  -c2(4)G(nt0KT -g]w  -  (x -g[2  -2C,«*)  -c2o»G(Ot0)( x  -g]  =  0 

(D.28) 

with  solutions 

w!(t0;^,x)  ■  41/CI(«^)G(Qt0)]{C1(^)  ♦  C2(0)G(ntoXx  -g  *  [Cfo)  +  2C3(^)G(QtoXx  -g]w}.  (D.29) 
Subtracting  Eq.  (D.28),  with  w  there  replaced  by  wl(t0;$,x),  from  Eq.  (D.15),  with  notation 

w(tc;  |x|,$,x)  *  t(t0\  |x|£,  |x|/c  +  x)  -  x  (D.30) 

in  the  latter  equation,  we  find,  using  Eq.  (D.25),  that  for  all  |x|  >  0 
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|[w(t0;  |x|,<»,T)  ♦  Q(t0,<M)]3  -  (w!(t0^,x)  ♦  0(to,$.x)]2|  <;  2(|x|/c)-1[l/C2((»G(Qt0)]rVG \Ot\ 


where 


(D.31) 


QM,x)  -  jqotf  ♦  c,«»G<ftt0Xx  -gyc^XXftg. 


(D.32) 


Hence  there  are  four  possibilities: 


Iimw(t0;  |x|,*,x)  -  wl(to;#,x), 
t»h*- 


(D.33) 


limw(t0;  |x|,$,x)  «  w:(t„;4,x), 
M-»- 


(D.34) 


lim  w(t0;  |x|,*,x)  «  4wl(t0;<M)  +2Q(t0,<M)). 


M-*- 


(D.35) 


and 


limw(t0;  |x|4,t)  «  4w;(t„;4,x) +2Q(t0,*,x)], 


(D.36) 


corresponding  id  the  four  roots  of  Eq.  (D.15).  Hie  correct  choice  is  easily  determined  by  noting  that 


x  ♦  lim  wI(to;$,x)  -  -x/C,(*)  -  x  ♦  limHw!(t0;4,x)+2Q(t0.*,x)]} 


(D.37) 
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to  that  these  three  choices  for  the  limit  yield  (for  x  >  0) 


Urn  Urn  I'd.;  |x|i,  |x|/c  ♦t)  -  t  «■  lim  w(t0;  |x|>,x)  -  -x/C.(6)  <  0 
«-»<r  Nh*-  i -»<r 


(D.38) 


yet,  by  Eq.  (5.16),  foe  iterated  limit  must  be  £  to  £  0.  On  foe  other  hand 


t  ♦  lim  wl(to;$,T)  -  t/C ,(♦); 


(D.39) 


hence  foe  correct  choice  is  Eq.  (D.33)  and  so 


lim  t  (t0;  |x|i,  |x|/c  ♦  x)  «  x  «■  wl(t0;*,x) 

Pr4- 


(D.40) 


as  claimed  in  Eq.  (D.12). 

(ii)  The  continuity  of  the  map  t„  ~  fl(tc;$,x)  at  to  =  0,  hence  on  [0,  f0"(x)] ,  follows  from 
Eqs.  (D.39)  and  (D.40),  to  wit. 


U®  *„(*„;♦•  t)  -  x  ♦  lim  w_4(tc;*,x)  -  x/C,(^)  -  t(0;6,x). 


(D.41) 


(iii)  Since  t  (t0;  |x|ft,  |x|/c  +  x)  -*  t^(tB;^,x)  to-pointwise  on  compact  [0,  T0"(x)]  and 
*o  **  Is  continuous  there,  then  by  Dini’s  Theorem  it  is  sufficient  to  show  that  foe  map 

1*1  ~  ^(t0;  lxl*-  l*(/c  +  x)  is  monotonic  for  each  tD  e  [0,  f  “(x)J;  in  fact,  we  show  that  there 
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exists  X  >  0  such  that  V,  ( |x|  >  X  ->  the  aforementioned  map  is  increasing.  To  that  end, 
we  apply  -  as  in  Theorem  C.4  -  the  Implicit  Function  theorem  to 

!Hi|x|,w)  ■  w  *  |x|/c  ♦  (fxl/cMsin2^  +  [|x|'*Z(to;  w  +x)  -  cos4JJ{“  (D.42) 

to  conclude  via  Eq.  (D.14)  that  |x|  ~  w(t0;  |x|,<>,x)  -  ftt.;  1*1*.  i*|/c  ♦  t)  -  t  is  C1 
at  £  £  (0, ■»)  whenever  £  and  t0  e  [0,  T"(x)]  satisfy 

0  it  {sin^  [£*‘Z(to;t/(t0,  £i,  £/c  ♦  x))  -  cos$l2}w  ■  R(t0,£)  (D.43) 

(where  the  4  and  x  dependence  of  R  have  been  supressed);  and  further,  for  such  £  and  t0,  that 

0w/a|x|Xto;£,*.x)  -  (l/c){l  -  R(t„,£)  ♦  [R(t0.£)]-‘[R.(t0,£)  -  cos*]R.(t0.£)} 

t|1+  (l/c)[R(to,£)]-'[R.(t0,£)  -  cos<>](dZ/dtXto;  t  (to;  £  i,  £  /c  *  x))} 

(D.44) 

where 

R.(tc,  £  )  ■  £  -‘Z(t0;  ?(t0,  £  *,  £/c  +  x))  (D.45) 

and  die  denominator  of  dw/9|x|  is  strictly  positive  whenever  the  condition  of  Eq.  (D.43)  obtains  since 
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(l/c)|[R.  -  coa^VR  |  '|dZ/dt  |  £  1/c  •  1  •  vc  «  \Jc  <  1. 


(D.46) 


We  now  claim  that  there  exists  X  >  0  such  that  V,  €  je  e  (X-)->  R(t0,  *)  >  0  so  that  the 
condition  of  Eq.  (D.43)  obtains  and  the  map  |x|  -» w(t0;|x|,#,x)  is  C1  on  (X,  ~)  and  Eq.  (D.44)  is 
valid  there;  indeed,  this  positivity  follows  immediately  from  the  definition  of  R  in  Eq.  (D.43)  by  using 
proposition  D2.  Further,  when  |xj  >  X  then  foe  numerator  of  dw/9|x|  is  positive,  as  follows.  We 
first  note  that  If  cos#  *  0  then  R,(t„,  |x|)cos#  -  0  <  1,  while  if  cos#  >  0  then,  again  by  Proposition 
D.2,  there  exists  X'>  0  such  that  [af;wl  \x\  >  X  ->  R.(to,|x|)cos#  <  1  and,  WLOG,  we  may 
take  X'~  X.  Thus  in  either  case  we  have,  since  #  >  0,  that  R*  >  R2  cos2#  so 
{R^  -  2R.cos#  +  1}W  >  1  -  R.cos#;  hence,  noting  from  Eqs.  (D.43)  and  (D.45)  that 

R  -  (Ri  -  2R.coe#  ♦  1}“  (PAT) 


we  then  have 

R  -  R-‘[R.  -  cos#JR.  -  rl[R*  -  (R.  -  cos#)RJ  -  [1  -  R.cos#]/R  <  1  (D.48) 

so  that  the  numerator  of  dw/9|x|  is  positive  for  all  tQ  e  (0,  T"(x)]  whenever  |x|  >  X.  Thus  if 
|x|  >  X  then  (dw/3)xjXt0;  |x|  ,#,x)  >  0  for  an  t„  e  [  0,  T“(x)]  so  the  map 
|t,  |x|/c  + 1)  is  increasing  for  such  t0. 

(iv)  This  is  straightforward  to  verify  using  Eq.  (D.12).  ■ 
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Proposition  D.4 :  Let  1  =  (sio^,  0,  cos4)  with  0  <  9  ^  x  >  0,  and  tc  e  [0,  T0  (x)J .  Then 


lim  V(t  ;  f  (t  ;  |x|*  |x|/c  ♦  x))  -  V(t#;  tl(to;fx))  (D.49) 

w— 

tg-unUbnnly  on  (0,  T  “(x)] . 

Proof.  That  the  limit  is  as  stated  follows  from  Proposition  D.3(i)  and  continuity  of  the  map 
t  ~  V(t*;  t)  at  t  *  t  (tD;|x|ft,  |x|/c  +  x),  which  continuity  follows  from  Eq.  (5.16)  since  x  >  0. 

Further,  since 

v(t0;iL(to ;♦.*))  -  v.o  -  2rKxaga'.(tB*,t)  -  U>  0>.so) 

then,  since  |G(Qto)|  £  1  for  all  to  c  [0.  T  "(x)],  we  have  for  all  such  tc 

|V(to;il(to;4,x))  -  V(to; f(to;  |x|t,  |x|/c  ♦  x))|  S  2T'l\l'jt0^,x)  -  |x|i.  |x|/c  «•  x)| 

(DJI) 


so  by  Proposition  D.3(iii)  the  convergence  in  Eq.  (D.49)  is  to>uniform  on  [0,  T“(x)].  ■ 


It  is  easy  to  compute,  using  Eqs.  (D.12)  and  (D.50),  that,  when  te  *  0, 
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(Va/c)[l  -  2r-‘G(£it0KT  -  g]. 


if  $  *  x/2 


V(t#;fl(t0;#,x)yc  - 

(l/co*)(l  -  C,(4){1  ♦  PCj(4VC1j(4)]G(QI0Xt  -  g}w).  if  o  <  ♦  <  * 

(D.S2) 

and  of  course 

V(0;ll(0;4,x)yc  -  wjc.  (D.53) 

Since 

[2CJ(tyC1a(4)lG(Ql0XT  -g  S  4(v,/c)r1G(Qto)7TG(Qto)]-,/ll  -(vycjcosf]  <  4{yjc)  <  1 

(D34) 

then  the  second  line  of  Eq.  (DJS2)  is  given,  correct  to  first  order  in  \Jc,  by 

y(g£.(t.))fc  -  (yjc)[i  -  2r,G(Qt0xx  -gj  (0  <  ♦  <  my,  (DJ5) 

in  fact,  this  last  expression  is  exactly  true  for  +  =  m,  as  Eq.  (D.52)  shows,  and  for  t0  =  0  as  well,  as 
Eq.  (D.53)  shows.  In  Section  V  and  beyond,  we  take  Eq.  (D.55)  (Eq.  (5.23))  to  be  our  expression  for 
V(g?l(t0;4,x)yc  in  an  cases. 
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Proof  of  Theorem  5.1 :  By  Proposition  D.2,  there  exists  X  >  0  such  that 

1*1  *  X-*  |x|*‘Z(t0;f(t0;  |x|i,  |x|/c  «•  t))  <  1.  Define  functions  k^.  and  h  as 

follows: 

00  tJ:[0,f,T(x)]  x  [*>)  ->  [-2/sin*,  2/sii*]  by 

1*1)  *  A.(t0;t'(t0;  |x|i,  jx|/c  ♦  t)) 

-  (1/sUm^) |x | ->Z<t0; I '(V.  |x|*  |x|/c  ♦  t))  -  cot* 
with|*,(t0;  |x|)|<  (l/sin$)  ♦  coty  <  2/sin^; 

(b)  fy  [0,  T  “(x)]  x  (*»  -*  [-v^c,  \Jc]  by 
M*«;  l*D  *  VOo'-^O,;  |*|t,  |x|/c  ♦  t))/c; 

(c)  for  p,  q  e  0,  p  i  0.  h^:  [0,  T0"(x)]  x  [-2/sin*,  2/sin*]  ->  *  by 

hM(to,k)  «  G'tfJtyCKQt^'d  ♦ 

wi*  iA„(t„**>i *  <  °°  * 

By  Eq.  (5.16),  the  maps  t„  ►*  Z  (t„;  ^)  and  t0  V(t0;  V)  art  continuous  at  tc,  hence  so  are 
*«(*.!  1*1).  1  *  1 A  further  since  G  and  G'  are  continuous  on  [0,  T“(t)J,  then  so  is  A  on  its  domain 
HnaDy,  by  Propositions  D.3  and  D.4  and  Eq.  (D55),  |x|)  converges  uniformly  to  -cot*  on 
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[0,  f0"(T)J  as  |x|  ->  •»  and  *^t0;  |x|)  converges  uniformly  to  (vyc)[l  -  2T'lG(iitJ(r  -  t^)]  (to 
first  order  in  sjc)  on  (0,  T"(x)]  as  jxj  ->  The  results  of  this  theorem  that  follow  by  applying 
Lemma  D.l  to  E,(|x|i,  |x|/c  ♦  x)  and  B^jxji,  |x|/c  ♦  x)  as  specified  by  Eqs.  (S.l)  -  (S.3).  ■ 
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